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These lectures contain an introduction to instantons, calorons and dyons of the Yang-Mills gauge theory. 
Since we are interested in the mechanism of confinement and of the deconfinement phase transition at some 
critical temperature, the Yang-Mills theory is formulated and studied at nonzero temperatures. We introduce 
"calorons with a nontrivial holonomy" that are generalizations of instantons and can be viewed as "made of" 
constituent dyons. The quantum weight with which these calorons contribute to the Yang-Mills partition function 
is considered, and the ensuing statistical mechanics of the ensemble of interacting dyons is discussed. We argue 
that a simple semiclassical picture based on dyons satisfies all known criteria of confinement and explains the 
confinement-deconfinement phase transition. This refers not only to the SU{N) gauge groups where dyons lead 
to the expected behaviour of the observables with A*', but also to the exceptional G(2) group whose group 
center, unlike SU{N), is trivial. Despite being centerless, the G{2) gauge group possesses confinement at low 
temperatures, and a l'^* order deconfinement transition, according to several latest lattice simulations, indicating 
that confinement-deconfinement is not related to the group center. Dyons, however, reproduce this behaviour. 
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1. WHAT IS 'VACUUM'? 

Quantum field theory deals with fields fiuctu- 
ating in space and time. If the fields are free, 
i.e. non-interacting, their fluctuations arc plane 
waves with any momenta, called zero-point os- 
cillations. The ground state of a system of free 
fields with no sources (the 'vacuum') is not a zero 
field but an infinite set of plane waves - just as 
the ground state of a one-dimensional quantum- 
mechanical oscillator corresponds not to a parti- 
cle lying at the bottom of the potential but to a 
particle distributed about the minima and having 
the energy The ground state of a free quan- 



tum field theory has the energy ^ ^ where one 
sums over all eigenfrequencies of the free fields; 
in a 3+1-dimensional theory this sum diverges as 
the fourth power of the cutoff momentum. 

However, if the fields are interacting, and the 
nonlincarity is strong enough such that it is not 
tractable by perturbation theory, their can be 
large field fluctuations in the vacuum, that can- 
not be reduced to weakly perturbed zero-point 
oscillations. Those large vacuum fluctuations de- 
termine the content of the quantum field theory, 
but in a strong interacting case it is difficult to 
describe them mathematically. One has either 
to rely on some approximate methods, or study 
the dominant fiuctuations numerically, or hope 
to gain some insight from similar, e.g. supersym- 
metric field theories where, because of additional 
symmetry, exact analytical methods can be devel- 
oped. Confinement is a very remarkable and still 
an unexplained phenomenon, therefore it may be 
helpful to use all known methods to understand 
its microscopic mechanism. 

There are two basic approaches to study the 
vacuum in quantum theory: one is the Hamilto- 
nian or Schrodinger approach, the other is the 
path integral or Feynman approach. We start 
from simple quantum mechanics of a particle in 
a 1-dimensional potential V{q) where q is the co- 
ordinate, to remind the connection between the 
two approaches. 

-2 

Let the Lagrangian he L = ^ V{q) and the 

energy H + V{q). To find the (quantized) 

energy levels En and the stationary wave function 
il>n{<l) one solves the Schrodinger equation: 

Another way to find all energy levels is to com- 
pute the partition function of a system, as func- 
tion of temperature T°. On the one hand, the 
partition function is defined as a sum of Boltz- 
mann exponents for all energy levels, 

Z = ^e~-^ . (2) 

n 

On the other hand the partition function can be 
presented, according to Feynman, as a path in- 
tegral over all particle's trajectories periodic in 



3 



time, with the period T = 



(3) 



r'?(T)=go 
= jdqo I Dq{t) 



exp — — / dt 







-V{q{t)) 



A generic path integral is illustrated in Fig. 1. 




Path integral over periodic trajectories 



minus sign in the exponent of Eq. (3). However, 
such a trajectory is unique, and it comes with a 
vanishing weight or entropy. If a trajectory os- 
cillates somewhat about the minimum such that 
both the potential and kinetic energies are not 
too large, one finds very many such trajecto- 
ries, therefore they have a much larger statistical 
weight. The outcome of this fight - between min- 
imizing the energy and maximizing the entropy 
- is well known: typical trajectories contributing 
most to the path integral are of the type shown 
in Fig. 2. This is why the ground-state energy is 
not zero but Eq — plus corrections from the 
deviation of the potential from the quadratic one. 



q(t) 




a 2a 3a ... 



T= IIT 



Figure 1. Feynman path integral in Quantum Me- 
chanics, i.e. in d = -I- 1 quantum field theory. 



Figure 2. A typical trajectory in a potential well 
exhibits zero-point oscillations giving rise to the 
vacuum energy Eq = 



One can visualize typical trajectories contributing 
to the path integral (3) by discretizing it and com- 
puting a big number of ordinary integrals over co- 
ordinates at discrete intermediate times tr, = na: 



N 

Z = M lim I I / dqn e 

n=l •' 



(4) 



*S'discr — ^ ^ ' 



"T- ( q{tn)-q{tn-l) 



-V{q{tn)) 



In order to cut out the vacuum state with the 
lowest energy Eq one has to take the limit of large 
observation time T — > oo, or small temperature 
T° ^ 0. At first sight it may seem that the lead- 
ing contribution to the partition function (3) will 
be from a trivial trajectory corresponding to a 
particle lying all the time T at the bottom of the 
potential well: it minimizes both the potential 
and the kinetic energies both coming with the 



For more complicated potentials, e.g. the double- 
well potential typical trajectories are more com- 
plex, see Fig. 3. 

If the barrier between the potential wells is high 
and broad, the particle will oscillate for a long 
time in one well but eventually will travel to the 
other one, where it will oscillate again and eventu- 
ally tunnel back to the first well, and so on. The 
tunneling parts of the trajectory are called, in 
modern language, instantons and anti-instantons. 
The definition is: Instantons are classical tunnel- 
ing trajectories q{t) with minimal action, satis- 
fying the equation of motion (here: the Newton 
law) 



5S 
W) 



= 



or mq = 



dV 



(5) 



As the particle transits from one well to another 
it still experiences zero-point oscillations which, 
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Figure 3. A typical trajectory in a double-well 
potential exhibits zero-point oscillations about 
tunneling trajectories - instantons and anti- 
instantons. 



however, are now not around zero but about a 
nontrivial background trajectory. The times of 
transition are random: in order to get physical 
quantities (like the level splitting in the double- 
well potential) one has to integrate over all 'posi- 
tions' of instantons and anti-instantons in time. 

To visualize instantons in computer simula- 
tions, one can start computing the partition func- 
tion (4) on a Id Euclidean lattice by means of, say, 
the Metropolis algorithm. Stop it when the sys- 
tem reaches equilibrium and draw the trajectory 
the computer is currently working with: it will 
be something like presented in Fig. 3. While one 
will not see directly a smooth trajectory interpo- 
lating between the two potential wells (because 
of the inevitable zero-point oscillations around it, 
producing Gaussian noise) , in Quantum Mechan- 
ics which is a (O-l-l)-dimensional quantum field 
theory it is not difficult to smear out the quan- 
tum noise and reveal instantons. In (3-1-1) dimen- 
sions smearing out the normal quantum fluctu- 
ations about smooth instanton configurations is 
not a straightforward task. A typical configura- 
tion of the gluon field in a lattice simulation of the 
(3-|-l)-dimensional Yang-Mills theory is shown in 
Fig. 4 borrowed from J. Negele et al. [1]. It takes 
some effort to 'cool' the configuration down to 
something smooth. Eventually one sees instan- 
tons and anti-instantons with a naked eye. 

Fig. 4, upper row, is a direct analogue of Fig. 3 




Figure 4. Instantons and anti-instantons in Ad 
Yang-Mills theory are revealed after zero-point 
oscillations are smeared out. Upper row: a typi- 
cal full configuration of the gluon field from lat- 
tice simulations [1] in the {z,t) hyperplane with 
(x, y) fixed. Lower row: the same configuration 
after smearing clearly shows that quantum noise 
of the upper row is around 3 instantons and 2 
anti-instantons. Left column: action density as 
function of Euclidean time and one spatial coor- 
dinate. Right column: topological charge density 
of the same configuration. 



but for a quantum field theory in (3+1) dimen- 
sions. 

2. INSTANTONS IN YANG MILLS 
THEORY 

2.1. Yang— Mills theory in Hamiltonian for- 
mulation 

A fundamental fact [2,3] is that the poten- 
tial energy of the gluon field is a periodic func- 
tion in one particular direction in the infinite- 
dimensional functional space; in all other direc- 
tions the potential energy is oscillator-like. This 
is illustrated in Fig. 5. The nontrivial form of 
the potential energy implies that instantons exist 
in the Yang-Mills (YM) theory, corresponding to 
tunneling through potential energy barriers. 
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Figure 5. Potential energy of the gluon field is 
periodic in one direction and oscillator-like in all 
other directions in functional space. 



To observe the periodicity, we take the Aq = 
gauge, called Weyl or Hamiltonian gauge. In 
Euclidean space it is more appropriate to denote 
the time component of the YM field by A4 , so we 
use the = gauge. For simplicity we start 
from the 5/7(2) gauge group, a = I, 2, 3. 

The spatial YM potentials (x, t) can be con- 
sidered as an infinite set of the coordinates of the 
system, where « = 1,2,3, a = 1,2,3 and x are 
"labels" denoting various coordinates. The YM 
action is 



S 



1 



(6) 



dt 



1 



d^x - 



1 



where Ef — A^ is the electric field strength 
and Bf = ieyfe {djA% - OkA'j + e''^''A]Al) is the 
magnetic field strength. 

Apparently, the first term in Eq. (6) is the 
kinetic energy of the system of coordinates 
{A°(x, t)} while the second term is minus the po- 
tential energy being just the magnetic energy of 
the field. Upon quantization the electric field is 
replaced by the variational derivative, Ef{x) 
—ig'^S/SAf{x), if one uses the 'coordinate' rep- 
resentation for the wave functional. The func- 
tional Schrodinger equation for the wave func- 



tional '5 [A" (a;)] takes the form 



n^[A, 



2 {6Afix)y 



^{Bnx))'\MA^]^£nA.] (7) 



where £ is the eigenenergy of the state in ques- 
tion. The YM vacuum is the ground state of the 
Hamiltonian (7), corresponding to the lowest en- 
ergy £. 

Not all solving formally the Schrodinger 

equation describe physical states, however, but 
only those that are gauge invariant: 



f7(x) e SU{2). 



Af = n^fAifi + in'fd^n, 



(8) 



Example: perturbative wave functional 



At small gauge coupling constant ^ one 
can neglect the commutator term in the magnetic 
field, i3f « eijkdjAl, the Schrodinger equation 
(7) becomes that for an infinite set of coupled 
harmonic oscillators, and can be solved exactly. 
The ground state wave functional is given by a 
Gaussian that can be best written in terms of the 
coordinates A^ (p) labeled by the 3- momenta p of 
the fields: 



exp 
exp 



d^p 

(27r)3 

PtP] 



(9) 



' //i3,,3.. gf(x)i?f(y) 



describing 2 • {N'^ — 1) gluon waves (for the SU (N) 
gauge group) with transverse polarizations and 
energy |p|. The vacuum energy of these zero- 
point oscillations is the sum of the harmonic os- 
cillators' energies ^ ^ or, more precisely. 



(27r)3 



(10) 



where V is the 3d volume. We see that the en- 
ergy of the zero-point oscillations is quartically 
divergent. 
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2.2. Topology 

Let us introduce an important quantity called 
the Pontryagin index or the four-dimensional 
topological charge of the YM fields: 



(11) 



pa _ _ rpa 
■^fiv — 2 1"""^ a/3' 



The integrand in eq. (11) happens to be a full 



derivative of the four- vector K/j,: 



(12) 



Therefore, assuming the fields Ai are decreasing 
rapidly enough at spatial infinity, one can rewrite 
the 4-dimensional topological charge (11) as 

Qt = jd*x{dtKi+diKi) = jdtj^j^^Ki. (13) 

Introducing the Chem-Simons number 

Ncs = J rf'x (14) 

= Jd'^e'^>'(AtdjAt + ^e'^>'<^AtA'^jAij 

we see from Eq. (13) that Qt can be rewritten 
as the difference of the Chern-Simons numbers 
characterizing the fields at t = ±oo: 



Qt = Ncs i+oo) - iVcs(-oo). 



(15) 



The Chern-Simons number of the field has an 
important property that it c;an change by inte- 
gers under large gauge transformations. Indeed, 
under a general time-independent gauge transfor- 
mation, 

Ai^U^AiU + iU^diU, Ai = A^^, (16) 
the Chern-Simons number transforms as follows: 

Ncs Ncs + Nw (17) 
i 



+ 



^ J d^x e'^^dj Tr (diUU^Ak). 



The last term is a full derivative and can be omit- 
ted if Ai decreases sufficiently fast at spatial in- 
finity. A^vi' is the winding number of the gauge 
transformation (16): 



Nw = 



-^1 
247r2 J 



247r2 

Tr [{U'^diU){U'^djU){U^dkU)\ . 



(18) 



The SU{2) unitary matrix U of the gauge 

transformation (16) can be viewed as a map- 
ping from the 3-dimensional space onto the 3- 
dimensional sphere of parameters S^. If at spa- 
tial infinity we wish to have the same matrix U 
independently of the way we approach the infin- 
ity (and this is what is usually assumed), then the 
spatial infinity is in fact one point, so the mapping 
is topologically equivalent to that from to S^. 
This mapping is known to be nontrivial, meaning 
that mappings with different winding numbers 
are irreducible by smooth transformations to one 
another. The winding number of the gauge trans- 
formation is, analytically, given by eq. (18). As 
it is common for topological characteristics, the 
integrand in (18) is in fact a full derivative. For 
example, if we take the matrix U (x) in a "hedge- 
hog" form, U = exp[i(r • T)/r P{r)], eq. (18) can 
be rewritten as 



N. 



w 



TT f^^ dr 

integer 



sin-^ P = - 

TT 



sin 2P 





(19) 



since P{r) both at zero and at infinity needs to 
be multiples of tt if we wish U (r) to be unambigu- 
ously defined at the origin and at the infinity. 

Let us return now to the potential energy of 
the YM fields. 



(20) 



One can imagine plotting the potential energy 
surfaces over the Hilbert space of the coordinates 
v4"(x). It will be some complicated mountain 
country. If the field happens to be a pure gauge, 
Ai = iWdiU, the potential energy at such points 
of the Hilbert space is naturally zero. Imag- 
ine that we move along the "generalized coordi- 
nate" being the Chern-Simons number (14), fix- 
ing all other coordinates whatever they are. Let 
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us take some point ^^(x) with the potential en- 
ergy V. If we move to another point which is 
a gauge transformation of Af{x) with a wind- 
ing number A''^^, its potential energy will be ex- 
actly the same as it is strictly gauge invariant. 
However the Chern-Simons "coordinate" of the 
new point will be shifted by an integer Nw from 
the original one. We arrive to the conclusion 
first pointed out by Faddeev [2] and Jackiw and 
Rebbi [3] in 1976, that the potential energy of 
the YM fields is periodic in the particular coor- 
dinate called the Chern-Simons number. This is 
illustrated in Fig. 5. 

2.3. YM instantons in simple terms 

In perturbation theory one deals with zero- 
point quantum-mechanical fluctuations of the 

YM fields near one of the minima, say, at A'^cs = 
0. The non-linearity of the YM theory is taken 
into account as a perturbation, and results in se- 
ries in where g is the gauge coupling. In such 
approach one is apparently missing a possibility 
for the system to tunnel to another minimum, say, 
at A^cs = 1- Tunneling is a typical nonperturba- 
tive effect in the coupling constant. 

Instanton is a large fluctuation of the gluon 
field in imaginary (or Euclidean) time corre- 
sponding to quantum tunneling from one mini- 
mum of the potential energy to the neighbour one. 
Mathematically, it was discovered by Bclavin, 
Polyakov, Schwarz and Tyupkin [4] ; the tunneling 
interpretation was given by Gribov ^. The name 
'instanton' has been introduced by 't Hooft [6] 
who studied many of the key properties of those 
fluctuations. Anti-instantons are similar fluctu- 
ations but tunneling in the opposite direction. 
Physically, one can think of instantons in two 
ways: on the one hand it is a tunneling process oc- 
curring in time, on the other hand it is a localized 
pseudoparticle in the Euclidean space. 

Following the WKB approximation, the tun- 
neling amplitude can be estimated as cxp(— S'), 
where S is the action along the classical trajec- 



^Gribov made this comment at a seminar given by 
A. Polyakov in 1976 at the Petersburg Nuclear Physics 
Institute, that I attended. Polyakov acknowledges Gri- 
bov's tunneling interpretation of instantons in his famous 
paper [5]. 



tory in imaginary time, leading from the mini- 
mum at iVcs = at t = — oo to that at A^cs = 1 
at t = +CO. 

According to eq. (15) the M topological charge 
of such trajectory is Qt = 1- To find the best 
tunneling trajectory having the largest amplitude 
one has thus to minimize the YM action (6) pro- 
vided the topological charge (11) is fixed to be 
unity. This can be done using the following trick. 
Consider the inequality 

o<jd'x {F^.-F^^y 

= J d^x (2F2 - 2FF^ = Sg'^S - din'^Qr, (21) 



hence the action is restricted from below: 



/,2 



S > ^Qt = 



/,2 ■ 



(22) 



Therefore, the minimal action for a trajectory 
with a unity topological charge is equal to Stt^/ g^, 
which is achieved if the trajectory satisfies the 
self-duality equation [4]: 



pa _ pa 



(23) 



Notice that any solution of Eq. (23) is simulta- 
neously a solution of the general YM equation of 
motion D'^F^^, = 0: that is because the "second 
pair" of the Maxwell equations, DI^F^^ = 0, is 
satisfied identically. 

Thus, the tunneling amplitude can be esti- 
mated as 

A ~ exp(— Action) = exp Jd^^xF^^^ 



= expl--^ 



27r 

exp ( 

a. 



(24) 



It is non-analytic in the gauge coupling constant 
and hence instantons arc missed in all orders of 
the perturbation theory. However, it is not a rea- 
son to ignore tunneling. For example, tunneling 
of electrons from one atom to another in a metal 
is also a nonperturbative effect but we would get 
nowhere in understanding metals had we ignored 
it. 
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2.4. Instanton configuration 

To solve Eq. (23) let us recall a few facts about 
the Lorentz group S0{3,1). Since we are talk- 
ing about the tunneling fields which can only de- 
velop in imaginary time, it means that we have 
to consider the fields in Euclidean space-time, so 
that the Lorentz group is just S'0(4) locally iso- 
morphic to SU{2) X SU{2). The gauge poten- 
tials Af^ belong to the (|, ^) representation of the 
SU{2) X SU{2) group, while the field strength F^^ 
belongs to the reducible (f,0) + (0,1) represen- 
tation. In other words it means that one linear 
combination of ^),^ transforms as a vector of the 
left SU{2), and another combination transforms 
as a vector of the right SU{2). These combina- 
tions are 



(25) 



where rj, fj are the so-called 't Hooft symbols, see 
below. We see therefore that a self-dual field 
strength is a vector of the left SU{2) while its 
right part is zero. Keeping this experience in 
mind we look for the solution of the self-dual 
equation in the form 



'lij — '^aij, 



l-t-$(x2) 



(26) 



Using the formulae for the 't Hooft 's rj symbols 
one can easily check that the YM action can be 
rewritten as 



'"'7 



2\dT 8 ^ 



(27) 



This can be recognized as the action of the 
double-well potential whose minima lie at $ = 
±1, and T plays the role of "time"; p is an arbi- 
trary scale. The trajectory which tunnels from 1 
at T = —oo to —1 at T = +00 is 



tanh 



(D> 



(28) 



and its action (27) is 5 = Stt^/^^, as needed. 
Substituting the solution (28) into (26) we get 



1 -I- tanh 



a;2(x^ + p^) 
The correspondent field strength is 
4p2 



pa 



(x2 + p2)2 



lua 



2f7^ 

laa 



(29) 



(30) 



192p4 



P 



2\4 ■ 



and satisfies the self-duality condition (23). 

The anti-instanton corresponding to tunnel- 
ing in the opposite direction, from Ncs = 1 to 
Ncs = , satisfies the anti-self-dual equation, 
with F —F] its concrete form is given by 
Eqs.(29, 30) with the replacement fj ^ rj. 




Figure 6. Action density of the YM instanton as 
function of z, t at fixed x = y = 0. Instanton is 
a 'ball' in Ad possessing 0(4) symmetry. 



Eqs.(29, 30) describe the field of the instan- 
ton in the singular Lorenz gauge; the singularity 
of Afj^ at = is a gauge artifact: the gauge- 
invariant field strength squared is smooth at the 
origin. Formulae for instantons are more compact 
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in the Lorenz gauge but an instructive homework 
is to pass on to the Hamiltonian A4 = gauge and 
to check exphcitly that instantons indeed corre- 
spond to tunnehng from Ncs = 1 to Ncs = . 

2.5. Instanton collective coordinates 

The instanton field, Eq. (29), depends on an 
arbitrary scale parameter p which we shall call 
the instanton size, while the action, being scale 
invariant, is independent of p. One can obviously 
shift the position of the instanton to an arbitrary 
4-point Zfj_ - the action will not change cither. Fi- 
nally, one can rotate the instanton field in colour 
space by constant unitary matrices U. For the 
SU{2) gauge group this rotation is characterized 
by 3 parameters, e.g. by Euler angles. For a gen- 
eral SU{N) group the number of parameters is 
_ 1 (the total number of the SU{N) genera- 
tors) minus (N — 2)^ (the number of generators 
which do not affect the left upper 2x2 corner 
where the standard SU{2) instanton (29) is re- 
siding), that is AN — 5. These degrees of freedom 
are called instanton orientation in colour space. 
All in all there are 

4 (centre) + 1 (size) + (4A'' — 5) (orientations) 

= AN (31) 

so-called collective coordinates describing the 
field of the instanton, of which the action is inde- 
pendent. 

It is convenient to introduce 2x2 matrices 



or as 



such that 



(32) 



2ir"r?r 

then the instanton field with arbitrary center z^, 
size p and colour orientation U in the SU{N) 
gauge group can be written as 



A — Aa+a 



-tp''U[a-{x-z)+ -{x-z)^]U^ 

{X - Zf\p^ + {X- Z)2] 



Tr(rt'') = 



AO. 



2p'0-'fr,^{x^z)^, 
{x — z)'^[p'^ -I- (a; — z)2] ' 

O''^ = Tr {Uh^Ua^), 0''''0'"' 



(34) 



This is the explicit expression for the AN- 
parameter instanton field in the SU{N) gauge 
theory, written down in the singular Lorenz 
gauge. 

2.6. One-instanton weight 

The contribution of a saddle point - one iso- 
lated instanton - to the partition function is 
called the one-instanton weight. We have already 
estimated the tunneling amplitude in Eq. (24) but 
it is not sufficient: the pre- factor is very impor- 
tant. To the 1-loop accuracy, it has been first 
computed by 't Hooft [6] for the SU{2) colour 
group, and generalized to arbitrary SU{N) by 
C. Bernard [7]. 

The general field can be decomposed as a sum 
of a classical field of an instanton Ajj^{x, ^) where 
^ is a set of AN collective coordinates character- 
izing a given instanton (see Eq. (33)), and of a 
presumably small quantum field a^(a;): 



(35) 



There is a subtlety in this decomposition due to 
the gauge freedom. The action is 



Action = d^xFl, 



(36) 



j d^xa^W^^a^ + Oia"^). 



(33) 



25 

Here the term linear in drops out because the 
instanton field satisfies the equation of motion. 
The quadratic form W^u being a second-order dif- 
ferential operator, has AN zero modes related to 
the fact that the action does not depend on AN 
collective coordinates. This brings in a divergence 
in the functional integral over the quantum field 
which, however, can and should be qualified 
as integrals over the collective coordinates: cen- 
ter, size and orientations. Formally the functional 
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integral over gives 
1 



(37) 



which must be i) normahzed (to the determinant 
of the free quadratic form, i.e. with no back- 
ground field), ii) regularized (for example by us- 
ing the Pauli-Villars method) , and Hi) accounted 
for the zero modes. Actually one has to compute 
a "quadrupole" combination. 



det'VKdet(W^o+/"^) 
detW^o det(M^ + Ai2) 



(38) 



where Wq is the quadratic form with no back- 
ground field and /x^ is the Pauli-Villars mass play- 
ing the role of the ultraviolet cutoff; the prime re- 
minds that the zero modes should be removed and 
treated separately. The resulting one-instanton 
contribution to the partition function (normal- 
ized to the free one) is [6,7]: 



■'1— inst 
Zp.T. 



do{p) 



exp 



dp I d^^'-'^Udoip), (39) 



27r 

a.s{p) 

2ti 
as{p) 



2N 



(PP) 



h- (40) 



The fact that there are all in all 4A^ integrations 
over the collective coordinates Zfj_, p,U reflects AN 
zero modes in the instanton background. The 
combination of the UV cutoff p and the coupling 
constant given at that cutoff as{p) is the YM 
scale parameter 



A = p cxp — 



2n 



UN a,{p) 



(41) 



The numerical coefficient C{N) depends implic- 
itly on the regularization scheme used. In the 
Pauli-Villars scheme exploited above 



C{N) = 



4.60 exp(-1.68A^) 
7r2(Ar- l)!(iV-2)!' 



(42) 



This is how the cutoff-independent combination 
A (41) comes into being. The mechanism is called 
the 'transmutation of dimensions'. Henceforth all 
dimensional quantities will be expressed through 



0.1 0.2 0.3 



0.4 0.5 



Figure 7. Distribution in instanton sizes com- 
puted in two-loop order [8]. The fall-off at large 
p cannot be trusted as the two-loop approxima- 
tion breaks down there. 



A, the only scale parameter of the massless YM 
theory, which is very much welcome. 
Notice that the integral over the instanton sizes 
in Eq. (39) diverges as a high power of p at large 
p: this is of course the consequence of asymp- 
totic freedom. It means that individual instan- 
tons tend to swell! This circumstance plagued 
the instanton calculus for many years. 



YANG MILLS THEORY 
NONZERO TEMPERATURES 



AT 



In the previous section we have introduced in- 
stantons at zero temperature, that is in the Eu- 
clidean space. We attempt, however, not only 
to understand confinement at zero temperature 
but also deconfinement at a high temperature. 
Therefore, we have to introduce the formalism for 
studying YM theory at a nonzero temperature [9] . 

We recall Feynman's representation for the par- 
tition function (see section 1): 



.4 



n 



q(P)=q / rl3 ' 

dq I Dq{t) cxp -/ di^Eudidfe, <?] 

9(0)=« V "'0 



In YM theory the role of coordinates q is played 
by the spatial components of the gluon field 
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DA,{t,x) 



Ai{0,x)=Ai(x) 

13 

dt I (fx 



Af 



However, in a gauge theory one sums not over all 
possible but only over physical states, i.e. satisfy- 
ing Gauss' law. In the absence of external sources 
it means that only those states need to be taken 
into account that are invariant under gauge trans- 
formations (see subsection 2.1): 



Ai{x) ^ [A{x)] 



Qix) 



n{x) 



Juia(x)t'' 



= n{x)^ A,{x) n{x)+in{x)^ ^^n{x). 

To restrict the summation to physical states, one 
projects to the physical i. e. gauge invariant states 
by averaging the initial and final configurations 
over gauge rotations. [This is as if in a quantum 
mechanics problem with a central-symmetric po- 
tential we would like to restrict the summation to 
spherically-symmetric states only.] We have 



phys states 



0n\ 



n) 



n ■' ■' 

/■ /■Ai(/3,cc)=Ai(cc)"2(x) 

= pni^2{x)DAi{x) DAi{t,x) 

J J Ai(0,x)=Ai(x)"i<^'''> 



exp< 



dt Id^i 



(Ai)+{Bty 



Renaming the initial field A^^ — > Aj and intro- 
ducing the relative gauge transformation fl{x) = 
^2(2;) ^l\{x) one can rewrite this as 



r fAdx)"^-^ 
^phys =pn{x)DA,{x) DA,{t, x) e 

J J Ai{x) 



-S[A,] 



.(43) 



A more customary form of the YM partition 
function is obtained if one introduces, instead 
of Q.{x), a new variable Ai(t,x), e.g. as = 
iexp(— ita;"t'')9texp(ifa;"f"). Then the partition 



function can be presented as a path integral over 
strictly periodic gauge fields: 

.Ephys = jDA,e^v\-^ jd''xF;^F;}^,{AA) 

A^{t,x)=A^{t + p,x), /?=^. 

An important variable is the Polyakov line or 

holonomy; it is the path-ordered exponent in the 
time direction, hence it can depend only on the 
space point x: 

L(x) = V exp dtAi{t,-x.)^ . 

In the first formulation (43) it is nothing but the 

group element n{x) by which the final Ai{x,P) 
can differ from the initial Ai(x,0): 

L{x) = Vl{x). 

Under space-dependent gauge transformations it 
transforms as i — > U~^LU. The eigenvalues of 

i(x) arc gauge invariant; for the SU{N) gauge 
group wc parameterize them as 



(45) 



L = diag (e 



(46) 



+ . . .+ j2N = 0, and assume that the phases of 
these eigenvalues are ordered: < < ■ ■ ■ < 
Mjv < = + We shall call the set of 

N phases {nm} the "holonomy" for short. Ap- 
parently, shifting ju's by integers does not change 
the eigenvalues, hence all quantities have to be 
periodic in all /u's with a period equal to unity. 

The holonomy is said to be "trivial" if L be- 
longs to one of the elements of the group cen- 
ter Zjsf. For example, in SU(3) the three trivial 
holonomies are 



Ml 



Ml = = = 
2 1 1 



Ml 



-3,M2 
1 

-3,M2 = 



3,M3 
1 

-3,M3 



L = l3, 

L = I3, 



L 



1.3. 



Trivial holonomy corresponds to equal yu's, mod- 
ulo unity. Out of all possible combinations of n's 
a distinguished role is played by equidistant fx's 
corresponding to Tr L = 0: 

1 1 



conf __ 



2N 



m 



(47) 
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For example, in SU (3) it is 

1 1 

Ml = ^3'^2 = 0,^3 = - 

=^ L — diag (^e^~ , 1, 



(48) 



We shall call it "most non-trivial" or "confining" 
holonomy. 

Immediately, an interesting question arises: 
Imagine we take the YM partition function (43) 
and integrate out all degrees of freedom except 
the eigenvalues {/im} of the Polyakov loop i(x) 
which, in addition, we take slowly varying in 
space. What is the effective action for /x's? What 
set of /z's is preferred dynamically by the YM sys- 
tem of fields? 



First of all, one can address this question in 
perturbation theory: the result for the potential 
energy as function of /i's is [9,10] 



ppc 



(2^)2T3 



N 



y^(Mm-A^n)^[l-(A^m-A^«)]^(49) 



It is proportional to (by dimensions) and has 
exactly N zero minima when all /i's are equal 
modulo unity, see Fig. 8. Hence, pp<=''* says that 
at least at high temperatures the system prefers 
one of the trivial holonomies corresponding to 
the Polyakov loop being one of the N elements of 
the center Zj^. However, terms with gradients of 
/i's in the effective action become negative near 
"trivial" holonomy, signalling its instability even 
in perturbation theory [11]. 
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Figure 8. The perturbative potential energy as 
function of the Polyakov line for the SU{2) (top) 
and SU{3) (bottom) groups. It has minima where 
the Polyakov loop is one of the N elements of 
the center Z^r and is maximal at the "confining" 
holonomy. 



Figure 9. Dyon-induced nonperturbative poten- 
tial energy as function of the Polyakov line for the 
SU{2) (top) and SU{3) (bottom) groups. Con- 
trary to the perturbative potential energy, it has 
a single and non-degenerate minimum at the con- 
fining holonomy corresponding to Tr L = 0. 
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It is interesting that in the supersymmetric 
Af = 1 version of the YM theory (where in ad- 
dition to gluons there are spin-i gluinos in the 
adjoint representation) the perturbative poten- 
tial energy (49) is absent in all orders owing to 
fcrmion-boson cancelation, but the nonperturba- 
tive potential energy is nonzero. Moreover, it is 
known exactly as function of ^'s [12]: it has a sin- 
gle minimum at precisely the "most non-trivial" 
or "confining" holonomy (47). The result can be 
traced to the semiclassical contribution of dyons, 
which turns out to be exact owing to supersym- 
metry. 

In the non-supersymmetric pure YM theory, 
the dyon-induced contribution cannot be com- 
puted exactly but only in the semiclassical ap- 
proximation (this is what the rest of the paper is 
about), and the perturbative contribution (49) is 
present, too. We shall show below that a semi- 
classical configuration - an ensemble of dyons 
with quantum fluctuations about it - generates 
a nonperturbative free energy shown in Fig. 9. It 
has the opposite behavior of the perturbative one, 
having the minimum at the equidistant (confin- 
ing) values of the /u's. There is a fight between 
the perturbative and nonperturbative contribu- 
tions to the free energy [13]. Since the perturba- 
tive contribution to the free energy is ~ with 
respect to the nonperturbative one, it certainly 
wins when temperatures arc high enough, and the 
system is then forced into one of the vacua thus 
breaking spontaneously the symmetry. At 
low temperatures the nonperturbative contribu- 
tion prevails forcing the system into the confining 
vacuum. At a critical Tc there is a confinement- 
deconfinement phase transition. It turns out to 
be of the second order for N = 2 but first order 
for N = 3 and higher, in agreement with lattice 
findings. 

For the centerless G(2) group, there is a 1^* 
order deconfinement transition (also in agreement 
with the lattice) although it is not associated with 
any symmetry. 

4. CONFINEMENT CRITERIA 

Confinement, as we understand it today and 
learn from lattice experiments with a pure glue 



theory, has in fact many facets, and all have to 
be explained. For example, in a general SU{N) 
group one can consider "quarks" in various ir- 
reducible representations. From the confinement 
viewpoint all representations are characterized by 
the phase it acquires under the gauge transforma- 
tion from the group center. The representation is 
said to have "A^'-ality = k" if the phase is 
Let us formulate mathematically the main con- 
finement requirements that need to be satisfied: 

• the average Polyakov line in any A^-ality 
nonzero representation of the SU{N) group 
is zero below Tc and nonzero above it 

• the potential energy of two static color 

sources (defined through the correlation 
function of two Polyakov lines) asymptot- 
ically rises linearly with the separation; the 
slope called the string tension depends only 
on the iV-ality of the sources 

• the average of the spatial Wilson loop de- 
cays exponentially with the area spanning 
the contour; at vanishing temperatures the 
spatial ("magnetic") string tension has to 
coincide with the "electric" one, for all rep- 
resentations 

• the mass gap: no massless gluons should be 
left in the spectrum; the massive glueballs 
must exhibit the Hagedorn spectrum, i.e. 
an exponentially rising density of states as 
function of mass 

• no Stefan-Boltzmann law for the free en- 
ergy, F ^ N'^T^; instead the free energy 
should be 0{N^) and very small below T^. 

5. SADDLE-POINT GLUON FIELDS 

The dynamics of the quantum YM system is 
governed by the YM partition function (43) or 
(44). A reasonable question to ask is whether 
there are field configurations A^(x, f) that give 
relatively large contributions to the partition 
function. Naturally, when we speak about certain 
field configurations we always imply that there 
are zero-point quantum fluctuations around, as 
one cannot stop them. 
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The privileged field configurations are appar- 
ently those that are local maxima of the weight 
for a given configuration, e^'^I'^l or local minima 
of the YM action S[A]. They are called saddle- 
point fields and satisfy the non-linear Maxwell 
equation, 



SS 







or d^^'f;:^ 



0, 



and have finite action S — J d'^x F^^^jF^^ < oo. 

5.1. Topological classification of classical 
configur at ions 

According to Gross, Pisarski and Yaffe [9], 
classical configurations with a finite action can 
be classified by the following numbers: 

1. Topological charge 

2. Holonomy, more precisely the gauge- invariant 
eigenvalues of the Polyakov loop at spatial infin- 
ity: 

eigenvalues of L — V exp i / dtA4{x,t) 

3. Magnetic charge, more precisely the gauge- 
invariant eigenvalues of the chromo- magnetic fiux 
at spatial infinity: 

/ 1 \ 
B = X -1 , etc. 

5.2. Generalization of standard instantons 

to T ^ 

In section 2 we have introduced the standard 
Belavin-Polyakov-Schwarz-Tyupkin (BPST) in- 
stanton that is a saddle-point YM field configura- 
tion at zero temperature. It is easy to generalize 
it to nonzero temperatures: to that end one has 
to make the fields periodic in the Euclidean time 
direction, with the period /3 = y. This has been 
done by Harrington and Shepard [14]. According 
to the classification of the previous subsection one 
obtains a configuration that has Qt ~ 1 and triv- 
ial holonomy, as the instanton field falls off fast 



enough at spatial infinity, see Eq. (34). For the 
same reason the magnetic charge of such configu- 
ration is zero. It has not 0(4) but only the 0(3) 
symmetry: it is a 'ball' in 3c? but squeezed by 
the periodic conditions in the time direction, see 
Fig. 10. 




Figure 10. A typical action density of the periodic 
instanton with trivial holonomy as function of z, t 
at fixed X = y = 0. It is 0(3) symmetric in space. 



The 1-loop weight of the periodic instanton 
(also called the caloron) has been computed by 
Gross, Pisarski and Yaffe [9] and has the form 
similar to Eq. (40): 



dz4 I d'^z 



d^U ^ 



dp C(2) / 877^ ^ ^ 



•(Ap)*(7rpr)t e-^^'^''^)' at upT > 1. 

Important: large-size calorons are suppressed by 
the Debye screening of the chromoelcctric field, 
represented by the last factor in Eq. (50). 

If the temperatures are high enough, the 
caloron sizes are, according to Eq. (50), strongly 
suppressed. Therefore, a 'dilute gas' of many 
calorons and anti-calorons is an approximate sad- 
dle point of the partition function, and one can 
easily evaluate it in that approximation. At lower 
temperatures the approximation fails, however 
one can still evaluate the partition function us- 
ing the variational principle developed in Ref. [8] 
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for zero-temperature instantons and applied to 
calorons in Ref. [15]. 

Do Harrington-Shepard calorons with trivial 
holonomy lead to confinement? In Fig. 11 we 
show the average < Tr L > as function of tem- 
perature. 



tential (49). The philosophy is then [13] that one 
has to admit that A4 is nonzero at low temper- 
atures, look for saddle-point configurations with 
a nontrivial holonomy, and check that they in- 
duce a nonperturbative potential that indeed has 
a minimum at a needed 'most nontrivial' value of 
^4. 



0.2 



300 



500 



Figure 11. The average Polyakov line (big 
dots) as function of temperature in MeV, for 
the vacuum built of periodic Harrington-Shepard 
calorons with trivial holonomy (D. Diakonov and 
V. Petrov, 1998, unpubHshed). The solid line 
is drawn approximately through lattice data for 
the SU{2) gauge group (the error bars are not 
shown) . 



It can be seen that it qualitatively follows the 
lattice data rising fast at T w A, however < Tr L > 
at low temperatures is close to but not strictly 
zero; there is no abrupt deconfincmcnt phase 
transition either. Therefore calorons with triv- 
ial holonomy do not and can not explain con- 
finement. The reason is that periodic instantons 
have A4 = at spatial infinity, hence a strict 
zero for < Tr L > needed for confinement can be 
achieved only through a miraculous cancelation, 
but it does not happen. 

The conclusion is that a modification of peri- 
odic instantons such that the field has A4 ~ ttT ^ 
at infinity, leading to Tr L = 0, would be helpful 
to obtain confinement. 

A conceptual difficulty, however, is that a 
nonzero A4 is forbidden by the perturbative po- 



6. BOGOMOLNY PRASAD SOMMER- 
FIELD (BPS) MONOPOLE 

BPS monopole or dyon is a static, time- 
independent self-dual solution of the YM equa- 
tion of motion Df^F^i, = having a finite action 
and topological charge, possessing both chromo- 
electric and chromomagnetic charges, and corre- 
sponding to nontrivial eigenvalues of the Polyakov 
line (or holonomy) at spatial infinity. 

These are gauge invariant characterizations of 
the solution, however to find it explicitly one 
needs to fix the gauge. We first use the so-called 
'hedgehog gauge' where the solution is 0(3) sym- 
metric and time-independent. Later we discuss 
BPS monopoles in other gauges. 

The key point is the nontrivial holonomy, i.e. 
nontrivial eigenvalues of the Polyakov line far 
away from the center of the BPS monopole. To 
ensure it, we choose the gauge where A4 is static 
and its modulus is constant at spatial infinity. 
BPS monopole is basically an SU{2) construc- 
tion, therefore we first build it for the SU{2) 
gauge group and then generalize it to higher-rank 
groups by embedding. 

We look for the monopole solution in the form 



Atix) 



A^x) 



E{r) 



1 - A{r) _ S(r) 



C{r) 



(51) 



where Ui is a unit vector in 3-space and A, B, C, E 
are functions of the distance r from the origin 
only. They are all dimensionless. 

There is a time-independent gauge transforma- 
tion which leaves the form of the above Ansatz 
unchanged. Indeed, if we perform a gauge trans- 
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formation 



the equations of motion: 



^ UA^U^+iUd^U\ U = exp {i{n ■ r)P(r)) , 

we remain inside the Ansatz but with the trans- 
formed functions: 



E" = E^, Dir) = 0, 



R" = R 



£2 + ij2 _ 1 



A Acos2P - Bsm2P, 

B Asm2P + B cos 2P, 

C C + 2rP', E^E. 

Let us introduce the electric (E) and magnetic 
(B) field strengths: 

pa _ C) Aa _ P) AO. ^ fO.bc Ab AC 

Et = F-, = DfAl Bf = \e,jkFfk. 

Substituting the potentials we get for the electric 
and magnetic fields: 



Ef 



-Da _ 



BE 
• r 

B' AC 



AE 



d E 



r r r r ar r 



+ riaUi 



r 

A^+B^-1 



J A' BC 

, +{Oai-nani)[ 

r I \rr 



The correspondent field strengths squared are 



d eV e^r-" 

-A 

dr r I r 



where we have parameterized the functions 
A, B, C as follows: 

A = i?cosa, B = Rsina, C = D + ra' . 



Since D{r) occurs to be zero one can choose the 
gauge transformation such that a{r) = 0, B{r) = 
0, C(r) = 0, and A{r) = R{r). A solution 
of these equations was first found (numerically) 
by Prasad and Sommerfield. Later Bogomolny 
showed that the solutions can be easily found an- 
alytically: 



A{r) = R{r) = 
E{r) 



vr 



1- 



sinh vr 1 
± (1 — urcothur) 

=F (vr-l + O (e 



0(e-'"'), r 



0, 
oo 



oo 



This analytical solution can be found in a simple 
way, by solving a first order equation. To see that, 
we notice that the energy functional satisfies the 
so-called Bogomolny inequality: 



/ d^r{Ef-Bff + -iz J d^rEfBf 
jd^r{Ef + Btf-%Jd\EtBt 



d^rEfB", 



The equality is achieved when the fields satisfy 
the Euclidean self-duality (or anti-self-duality) 
equations. 




pa 



pa 



or Ef = ±Bf. 



Under this gauge transformation a{r) — > a(r) + 
2P(r) while the functions R{r) and D(r) rest in- 
variant, hence only they appear in the gauge in- 
variant expressions. The Euclidean energy func- 
tional is 



{Etf + {Btf 



Varying this functional with respect to the 
gauge invariant functions D,E and R one finds 



It is easy to chock that the two solutions satisfy 
exactly these equations, so that the upper sign 
is in correspondence with the upper sign in the 
previous equation and vice versa. In both cases 
the energy of the solution can be computed as 

g -1 

9^ 



I 



dr 



d^rEfBf 
d {1-A^)E 



dr 



v>0. 
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Let us write down explicitly the electric and 
magnetic field strengths: 



Bf 



E? 



sinh vr 



- — coth vr 
r 



+nani 
±B?. 



1 



sinh^ vr 



We see that both electric and magnetic field 
strengths fall off at spatial infinity as 1/r^; this 
behaviour is characteristic of the fields of elec- 
tric and magnetic charges. That is why the BPS 
monopole is often called 'dyon'. 

6.1. Four dyons of the SU{2) gauge group 

(M, M,L,L) 
In SU{2) there are four dyons with all four pos- 
sible signs of the electric and magnetic charges. 
The two usual BPS dyons in the regular ('hedge- 
hog') gauge have the form: 



Al = TnaV^ivr), 
1 

^{z) = coth^; 



z 



At = 
R{z) = 



Z 

sinh 2; 



1 - R{vr) 



O(0e-^). 



The upper sign in A4 corresponds to the self-dual 
{E^ = Bf) and the lower sign to the anti-self-dual 
{Ef = -Bf) solution. We shall call them M- and 
M-monopoles, respectively. The magnetic field 
strength in the hedgehog gauge is given by two 
structures: 



-Bf = {8ai - UaUi) Fi{r) + UaUi i^2(r) 



where 

F2(r) 



-^R{vr) 
r dr 



R{vr)^{vr) 



(52) 



fl-coth(t;r)) = v^O (e— ) , 
inh(ur) \vr J 

R^ivr)-! 



— —v^(vr) = 
dr 



(53) 



sinh(t;r) 



l = -lWO(e--). 



If there is more than one monopole in the vac- 
uum it is impossible to add them up in the hedge- 
hog gauge: one has to "gauge-comb" them to a 
gauge where A^ has the same asymptotic value 
at spatial infinity for all monopoles involved, - 
say. along the third colour axis. It is achieved 
with the help of two unitary matrices dependent 
on the spherical angles 0, <j): 



S+{ 



s. 



■ d) 3 ■ 2 -4) 3 
— 7— T 'I — T 'I —T 



S+(n ■ t)S\_ = T'\ 

3 i^L^T^ J* 3 

e 2^ e 2 ^ , 



S-{n-T)S'. 



We shall gauge-transform the M-monopole field 
with S- and the M-monopole with S+. As the 
result their A4 components become equal: 



' = v^{vr)- 



1 



O e- 



On the contrary, the spatial components differ in 
sign. We write them in spherical components: 



0, 



Ae = ^{r' sin ct> + T^ cos <!>), 



A. 



^p-{T COS (j) — T sin( 



+ ^tan|T3. 



The azimuthal component of the gauge field has 
a singularity along the negative z axis, therefore 
we shall call it the "stringy gauge". The field 
strength, however, has no singularities. The elec- 
tric field both of M and M monopoles in the 
stringy gauge is 



E. 



M,M 



Ef — 
E$ = 

Es = 



2 ^ r2 2 

^(-T^ cos.^ + t2 sin,^) 
(e-"'') 

^(r^ sin<^ + T2 cos<^) 
(e-""-) 
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while the magnetic field is -Bj = ±Ei. We see 
that the 9, (j) components are non-Abelian but 
exist only inside the dyon core, whereas the ra- 
dial component is Coulomb-like and Abelian as 
it has only the 3*^ colour component. Therefore, 
M monopolc has Abelian (electric, magnetic) 
charges (++) whereas the M one has (H — ). 

There is a second pair of dyons: a self-dual one 

with the charges ( ) which we shall name L- 

monopole, and an anti-self-dual one with charges 
( — h) which we shall name L monopole. They are 
obtained from the above equations by replacing 
V — > 2'kT — v. One first transforms them from the 
hedgehog to the stringy gauge with the help of 
the unitary matrices 5+ and S- , respectively. As 
a result, they get the same asymptotics Ai{<xi) = 
(— 27rr -|- u -|- ^) To put the asymptotics in 
the same form as for M, M-monopoles one makes 
an additional gauge transformation with the help 
of the time-dependent matrix 

U = exp {-iirTx'^T^) . 

This gives the following fields of L, L monopoles 
in the stringy gauge: 

--3 

A^'^ = [{2ttT -v)^{\2ttT -v\r) -2ttT] — 



1\ 



E. 



L,L 



J'2(r) 3 
2 ' 



1 ll 
7^ 2 



2 



■T COS 



(/> + t2 sin 4>)U^ 



■ ^ i7(r 1 sin </> + cos (j)) W , 



B. 



= ±E, 



The 'profile' functions F12 are given by Eqs.(52, 
53), with the replacement v 2'kT—v. Wc> notice 
that "the interior" of the L,L dyons, represented 
by the 6, cp field components are time dependent. 
This is why in the true 3d case these objects do 
not exist! 

Properties of the four fundamental dyons of the 
SU (2) group are summarized in the Table: 
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e-charge 
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m-charge 
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action, -^-^ 
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V 


2'kT-v 


2-kT-v 


top-charge 
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27rT 
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2vT 


1 — 

27rT 


— 1 

2-kT ^ 



6.2. Dyons of higher-rank gauge groups 

Dyons or BPS monopoles are basically SU{2) 
constructions, therefore if one wants to construct 
dyon solutions for higher-rank gauge groups one 
has to embed the SU{2) construction of the pre- 
vious subsection into a higher-rank group. This, 
of course, can be done by infinitely many ways; 
the problem is to construct a set of 'fundamental' 
dyons from which all other types can be built as 
'bound states'. The problem has been addressed 
in Ref. [17] for all Lie groups. Here we shall con- 
centrate on the SU{N) series, and later add some 
remarks on the exceptional G{2) group which is 
interesting both from the mathematical and phys- 
ical points of view since, contrary to SU{N), this 
group has only a trivial center. 

In SU{N) there are exactly TV kinds of funda- 
mental self-dual dyons (and anti-self-dual anti- 
dyons) which can be described as follows. Let us 
introduce a basis of Cartan generators given by 
diagonal matrices 



C„ 



= diag(0,...,0,l,-l,0,...,0), 
1,...,N-1, 



(54) 



m ■■ 



where +1 is in the m place, and supplement this 
set by the matrix Cat = diag(— 1,0, ...,0, 1), such 
that X]m=i ^™ ~ 0- ^ types of dyon solu- 
tions are labeled by the holononiy or the set of 
eigenphases /j 's of the Polyakov line at spatial in- 
finity. In the gauge where A4 is time-independent 
it means the asymptotic field 

^4(00) = 27rTdiag(/xi,/i2, • . • ,Miv), (55) 
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N 

Ml < M2 < • • • Mat < Ml + 1) Mm = 0- 

m— 1 

We adopt cyclic-symmetric notations such that 
/Ltjv+i = /zi + 1 and introduce the differences 

J^m = Mm+1 - Mm, '^"^ ^ 

m— 1 

The dyon of the m*'^ kind is an object of the 
SU{2) whose generators are Cm playing the role 
of the Pauli matrix , and the two corresponding 
off-diagonal matrices playing the role of the Pauli 
matrices. At the origin its component is 

4"^(0) = 27rrdiag(/xi,... 

■■■l^'m—li ^ J ^ )Mm+2)"- I ) 

and its core (where the fields are strong and non- 

Abelian) has the size ^ l/(27rTi/m)- Beyond the 
core the fields are Abelian as they have only a 
diagonal component, and are Coulomb-like: 

±B(-)=E('") '"^ Cm^. (57) 

The upper sign for the magnetic field is for the 
self-dual dyons, and the lower sign is for the anti- 
self-dual anti-dyons. In the gauge where A4 is 
static, the last N^^ dyon is obtained through a 
time-dependent gauge transformation, like; the 'L' 
dyon of the previous subsection. It is sometimes 
called the Kaluza-Klein monopole; in a true 3d 
theory it does not exist. 

The action density is, however, time- 
independent for all N monopoles. The full ac- 
tion and the topological charge are temperature- 
independent: 

= —Vm, Q'"^ = Vm. (58) 

The full action of all N kinds of well-separated 
dyons together are that of one standard instan- 
ton: S'inst = 2-1^ jag. We also note that the to- 
tal electric and magnetic charge of a set of N 
dyons is zero, and that the topological charge 



is unity. Therefore, a set of N kinds of well- 
separated dyons have the quantum numbers (and 
the action) of one instanton. However, contrary 
to the standard instanton or its periodic general- 
ization, the holonomy is not trivial but arbitrary. 
Remarkably, there exists an exact classical solu- 
tion that preserves these properties even as one 
moves the N dyons close to each other, see the 
next section. 

For the G{2) group, there are three kinds of 
fundamental dyons (as in S'[/(3)) but a neutral 
combination ('the instanton') is obtained from 
four dyons: one of the three kinds has to be taken 
twice. 

In the semiclassical approach, one has first of 
all to find the statistical weight with which a 
given classical configuration enters the partition 
function. It is given by exp(— Action), times the 
(determinant)"^/^ from small quantum oscilla- 
tions about the saddle point. For an isolated 
dyon as a saddle-point configuration, this factor 
diverges linearly in the infrared region owing to 
the slow Coulomb decrease of the dyon field. It 
means that isolated dyons are not acceptable as 
saddle points: they have zero weight, despite 
finite classical action. However, one may look 
for classical solutions that are superpositions of 
N fundamental dyons, with zero net magnetic 
and electric charges. The small-oscillation deter- 
minant must be infrared-finite for such classical 
solutions, if they exist. 



7. INSTANTONS WITH NONTRIVIAL 

HOLONOMY 

The needed classical solution has been found a 
decade ago by Kraan and van Baal [18] and inde- 
pendently by Lee and Lu [19], see also [20]. We 
shall call them for short the "KvBLL instantons" ; 
an alternative name is "calorons with nontrivial 
holonomy" . The solution was first found for the 
SU (2) group but soon generalized to an arbitrary 
SU{N) [21], see [22] for a review. 

The general solution A^^^^^ depends on Eu- 
clidean time t and space x and is parameter- 
ized by 3A/' positions of N kinds of 'constituent' 
dyons in space xi, . . . ,xjv and their U{1) phases 
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Figure 12. Action density inside the SU{3) 
KvBLL instanton as function of time and one 
space coordinate, for large (top), intermediate 
{middle) and small (bottom) separations between 
the three constituent dyons. 



ipi, . . . ,tpN. All in all, there are AN collective co- 
ordinates characterizing the solution (called the 
moduli space), of which the action S'inst = ^n/as 
is in fact independent, as it should be for a gen- 
eral solution with a unity topological charge. The 
solution also depends explicitly on temperature T 
and on the holonomy fii, . . . , fj.^: 

^KvBLL ^ a; (t, x; XI, . . . , xjv, V'l, • • • , ^n; 

T,^l, . . . ,^Ar) . (59) 

The solution is a relatively simple expression 
given by elementary functions. If the holonomy 
is trivial (all /i's are equal modulo unity) the ex- 
pression takes the form of the strictly periodic 



0(3) symmetric caloron [14] reducing further to 
the standard 0(4) symmetric BPST instanton [4] 
in the T limit. At small temperatures but ar- 
bitrary holonomy, the KvBLL instanton also has 
only a small 0{T) difference with the standard 
instanton. 

One can plot the action density of the KvBLL 
instanton in various corners of the parameter 
(moduh) space, see Fig. 12. 

When all dyons are far apart one observes A'' 
static {i.e. time-independent) objects, the iso- 
lated dyons. As they merge, the configuration 
is not static anymore, it becomes a process in 
time. The reason why time dependence creeps 
in is that the field of one of the N kinds of dyons 
is time-dependent from the start (in the nota- 
tions of subsection 6.1 it is the L dyon). The 
time dependence of that particular dyon is pure 
gauge rotation for isolated objects, therefore the 
action density and other gauge-invariant quanti- 
ties are time-independent. However, as we merge 
together dyons and wish to support the property 
that the configuration remains to be a solution of 
the equation of motion, time dependence ceases 
to be a gauge rotation but enters the field in an 
essential way, see the next subsection. 

In the limiting case of a complete merger, the 
configuration becomes a Ad lump resembling the 
standard instanton. The full (integrated) action 
is exactly the same S'inst = ^n/as for any choice of 
dyon separations. It means that classically dyons 
do not interact. However, they do experience a 
peculiar interaction at the quantum level, which 
we discuss in Section 8. 

7.1. Explicit form of the KvBLL instanton 

for the SU{2) gauge group 

We follow in this subsection Ref. [18]. We con- 
sider a classical solution whose asymptotic A4 
component is A4 — > v^. In the notations of the 
previous subsection v = 47rT/ii. We introduce 
also a notation v = 27tT° — v. 

The KvBLL instanton is constructed as a com- 
bination of L and M monopoles, see subsection 
6.1. Let the L dyon be centered at a 3d point xi 
and the M dyon is at a point X2. The dyon sep- 
aration is ri2 = [xi — X2[ = irTp^. Here p is the 
size of the standard instanton (see subsection 2.4) 
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to which the KvBLL instanton reduces at v ^ 
or at T ^ 0. 

We introduce the distances from the 'obser- 
vation point' X to the dyon centers, r = x — 
Xi, s = X — X2. Correspondingly, r = |r|, s = 
|s|. We choose the separation between dyons to 
be along the third spatial direction, ri2 = e3ri2. 

The KvBLL instanton field is 



(60) 



T 1 _ 

+ $ Re [(?7i. (r'+ir^) (9, +iv,5,,4)x] , 

where r" are Pauli matrices, fj'^^ is the 't Hooft 
symbol (see subsection 2.4). The functions used 
are 

i\ I „i, „!, I shsh. 



i/, = - cos(27rTx"*) + chch 



2rs 



■0 = ^ H — 3^shsh H shch H shch, 

rs s r 

■0 \ s r 



t 



We have introduced short-hand notations for hy- 
perbolic functions: 

sh = sinh(sv), ch = cosh(sv), 
sh = sinh(rv), ch = cosh(rv) . 

The first term corresponds to a constant A4 com- 
ponent at spatial infinity (A4 « «v^) and gives 
rise to the non-trivial holonomy. One can see that 
is periodic in time with period l/T. A useful 
formula for the field strength squared is 

Tr F^,^^^, = 92^2 iQg^^ 

In the situation when the separation between 
dyons ryi is large compared to both their core 
sizes i (M) and = (L), the caloron field can be ap- 
proximated by the sum of individual BPS dyons. 
To demonstrate it, we give below the field inside 
the cores and far away from both cores. 

7.2. Inside dyon cores 

In the vicinity of the L dyon center Xi and far 
away from the M dyon (sv 3> 1) the field be- 
comes that of the L dyon. It is instructive to 
write it in spherical coordinates centered at Xi. 




Figure 13. The action density of the KvBLL 
caloron as function of z^x at fixed t = y = 0. 
At large separations r\2 the caloron is a super- 
position of two BPS dyon solutions (left: r\2 = 
1.5/T). At small separations they merge {right: 
ri2 = 0.6/r). 



In the 'stringy' gauge in which the A4 compo- 
nent is constant and diagonal at spatial infinity, 
the L dyon field is 



= — { - + 2tt - vcoth(vr) 
2 \ r 



Ar 



0, 



A^ - V- 



sin(27rTa;^ 



+ cos(27ra;^ — 



A 



L _cos(27ra; 



2 sinh(vr) 
+ sin(27rTa;4 



2 sinh(vr) 



3tan(6'/2) 



Here Ag^ for example, is the projection of A onto 
the direction xiq = (cos cos 0, cos sin 0, — sin 0). 
The <j) component has a string singularity along 
the z axis going in the positive direction. No- 
tice that inside the core region (vr < 1) the field 
is time-dependent, although the action density is 
static. At large distances from the L dyon cen- 
ter, i.e. far outside the core one neglects expo- 
nentially small terms O (e"^'') and the surviving 
components are 



1 



y 



tan? 



corresponding to the radial electric and magnetic 
field components of the L dyon (see subsection 
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6.1): 



Et 



B. 



L r—^oo 



1 



Similarly, in the vicinity of the M dyon and far 
away from the L dyon (rv :» 1) the field becomes 
that of the M dyon, which we write in spherical 
coordinates centered at X2: 



^4 



V coth(vs) 



sm0 r 



COSfl 



= V 



2i sinh(vs) 
coscf) + sincj) 



2 sinh(vs) 
whose asymptotics is 



+ r 



= 0, 



;tan(^/2) 
' 2^ 



iM 



tan^ 



1 

s 

Y 



T 

y 



^5 = a 



M 



We see that in both cases the L,M fields become 
Abelian at large distances, corresponding to (elec- 
tric, magnetic) charges (— , — ) and (+, +), respec- 
tively. The corrections to the fields are hence of 
the order of 1/ ri2 arising from the presence of the 
other dyon. 

7.3. Far away from dyon cores 

Far away from both dyon cores (rv ^1, sv ^ 
1; note that it does not necessarily imply large 
separations - the dyons may even be overlapping) 
one can neglect both types of exponentially small 
terms, O (e^''^) and O (e^'^^). With exponential 
accuracy the function x is zero, and the KvBLL 
field becomes Abelian: 



T3 

2 



where is 



ri2 



S3 . 



r + s - ri2 



r3 



if ri2 = ri2e3. 



It is interesting that, despite being Abelian, the 
asymptotic field retains its self-duality. This is 



because the 3'^ colour component of the electric 
field is 

Ef = diAl = didaln 
while the magnetic field is 

Bf = e,,kdjAl = dids In $as - Sizd"" In 

where the last term is zero, except on the line 
connecting the dyon centers where it is singular; 
however, this singularity is an artifact of the expo- 
nential approximation used. Explicit evaluation 
gives the following nonzero components of the 
field far away from both dyon centers: 



1 1' 

, vH 

2 V r s. 



(61) 



T3 / 1 ^ 1 \ / {ri2~r+s){ri2+r-s) 
2 \r s J \l iri2+r+s){r+s-ri2) 

In particular, far away from both dyons, A4 is 
the Coulomb field of two opposite charges. The 
azimuthal A^ component is that of the magnetic 
field created by a current going from Xi to X2. 

8. QUANTUM WEIGHT OF A NEU- 
TRAL CLUSTER OF N DYONS 

Remarkably, the small-oscillation determinant 
about the KvBLL instanton can be computed ex- 
actly; this has been first done for the SU (2) group 
in Ref. [23] and later generalized to SU{N) in 
Ref. [24]. The quantum weight of the KvBLL 
instanton can be schematically written as an in- 
tegral over 3A^ coordinates of dyons (the weight 
does not depend on the U{1) angles tpm, hence 
they can be integrated out): 



f , , n /47rAt4 

J dyii...dxN \/detg ( ^ jT 

^-P(-^) (Det(-A))-^_ 



N 



(62) 



where g is the full 4iV x AN metric tensor of the 
moduli space, defined as the zero modes overlap 
matrix, and Dot (—A) is the functional determi- 
nant over nonzero modes, normalized to the free 
one and regularized by the background Pauli- 
Villars method; ^ is the Pauli-Villars ultraviolet 
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cutoff and g"^ is the bare coupling constant de- 
fined at that cutoff. The Jacobian det g turns out 

to be a square of the determinant oi a,n N x N 
matrix G(i) such that y/detg = det G^^^ where 



1 



(63) 



+ 



1 



m-l| 



is a matrix whose entries are Coulomb interac- 
tions between dyons that are nearest neighbours 
in kind. The Coulomb interactions in the zero 
mode overlap matrix arise naturally from the 
Coulomb asymptotics of the dyon field (61), so 
it is quite simple to check that Eq. (63) is cor- 
rect at large separations. A nontrivial fact is 
that Eq. (63) is actually exact for all separa- 
tions between dyons, including the case when 
they strongly overlap like in Fig. 12, bottom. This 
has been first conjectured by Lee, Weinberg and 
Yi [25] and then proved to be indeed exact at all 
separations by a direct calculation by Kraan [26] 
and later checked in Ref. [27]. In the last paper 
it has been also shown that in the limit of trivial 
holonomy (fim = 0) or vanishing temperature the 
measure given by Eqs.(62, 63) reduces to the stan- 
dard instanton measure (39) written in terms of 
the conventional "center-size-orientation", which 
is a rather nontrivial but gratifying statement. 

The functional determinant over nonzero 
modes Det (— A) together with the classical ac- 
tion and the Pauli-Villars cutoff combine into the 
renormalized scale parameter Apy, times a func- 
tion of dyon separations, Apy and T [23,24]. It 
is a complicated function which, for the time be- 
ing, we approximate by its most essential part: a 
constant equal to exp (— PP'^''*), where PP'"'' is the 
perturbative gluon loop (49) in the background of 
a constant field A4 (55). This part is necessarily 
present in Det^^(— A) as most of the 3d space 
outside the instanton's core is just a constant A4 
background, and indeed the calculation [23,24] 
exhibits this piece which is the only one propor- 
tional to the 3-volume. 

Therefore, we write the weight of the KvBLL 
instanton i.e. a neutral cluster of N different-kind 



dyons as 

Wi « jdxi...dxNdetGif^ 

• exp(-PP-*(/xi,...,Miv)) 

where the fugacity / is 



(64) 



(65) 



The bare coupling constant g^ in the pre- 
exponent is renormalized and starts to "run" only 
at the 2-loop level not considered here. Eventu- 
ally, its argument will be the largest scale in the 
vacuum, be it the temperature or the equilibrium 
density of dyons. 

9. QUANTUM WEIGHT OF MANY 
DYONS 

In the vacuum problem, one needs to use not 
one but 0{V) number of KvBLL instantons as 
the saddle point. Solutions with the topological 
charge greater than 2 will be hardly ever known 
explicitly as their construction runs into the prob- 
lem of resolving the nonlinear Atiyah-Drinfeld- 
Hitchin-Manin-Nahm constraints. At present 2- 
instanton solutions characterized by a nontrivial 
holonomy have been found [28] but it is insuffi- 
cient. Nevertheless, the moduli space measure of 
an arbitrary number K of KvBLL instantons can 
be constructed despite the lack of explicit solu- 
tions, at least in the approximation which seems 
to be relevant for the large- volume thermodynam- 
ics, if not exactly. 

When one takes a configuration of K instantons 
each made of N different-kind dyons one encoun- 
ters also same-kind dyons for which the metric 
(63) is inapplicable. However, the case of identi- 
cal dyons has been considered separately by Gib- 
bons and Manton [29]. The integration measure 
for K identical dyons following from that work is, 
in our notations, 



W 



ident 



/dxi...rfx^detG^^^"J,, (66) 



G 



ident 



47rf„ - Efe^i T\^- 
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where the identity factorial is inserted to avoid 
counting same configurations more than once. 

As in the case of different-kind dyons, this re- 
sult for the metric can be easily obtained at large 
separations from considering the asymptotics of 
the zero modes' overlap. The four zero modes 
^^T'^l't = 1,2,3,4) of individual dyons arc given 



by the components of the field strength: 



FuK- The zero modes for the m*^ kind of dyon are 



normalized to its action, J Tr 
(see Eq. (58)) and hence depend on the holon- 
omy. Since the field strengths decay as the 
overlaps between zero modes arc Coulomb-likc, 
and only those that are nearest neighbors in m 
do interact. In fact, the diagonal components of 
the metric tensor also acqiiire Coulomb-like cor- 
rections since the action of individual dyons is ac- 
tually normalized to its asymptotic field A4 that 
gets Coulomb corrections from other dyons. 

However, in contrast to the different-kind 
dyons, it is not possible to prove that this ex- 
pression is correct at all separations. Moreover, 
such an extension of Eq. (66) is probably wrong. 
The metric for two same-kind dyons has been 
found exactly at all separations by Atiyah and 
Hitchin [30] by an undirect method: it is more 
complicated than what follows from Eq. (66) at 
K — 2 but differs from it by terms that are ex- 
ponentially small at large separations [31]. We 
shall neglect the difference and use the Gibbons- 
Manton metric at face value. The point is, 
Eq. (66) imposes very strong repulsion between 
same-kind dyons (as does the Atiyah-Hitchin 
metric), hence the range of the moduli space 
where the two metrics differ is, statistically, not 
frequently visited by dyons. We do not have a 
proof that all thermodynamic quantities will be 
computed correctly with this simplification: prov- 
ing or disproving it is an interesting and impor- 
tant problem. To remain on the safe side, one has 
to admit today that the metric (66) is applicable 
if the dyon ensemble is sufficiently dilute, that is 
at high temperatures. Nevertheless, physical ob- 
servables we compute have a smooth limit even 
at T — > 0. Therefore, it may well prove to be 
a correct computation at any temperatures, but 
this remains to be seen. 



It is possible to combine the metric tensors for 
different-kind (63) and same-kind (66) dyons into 
one metric appropriate for the moduli space of 
Ki dyons of kind 1, K2 dyons of kind 2,..., Kn 
dyons of kind N. It is a matrix whose dimension 
is the total number of dyons, that is a {Ki -|- . . . -|- 
Kn) X (i^i + ... + Kn) matrix: 



G„ 



47ri/„ 



E 



■m+l,fe I 



E 



5 



in .ri — 1 



k 



mi ^m—l,k\ 



Tlx, 
+ 2 



■rm ^m-\-l,j 



T|x„ 



(67) 



where x^i is the coordinate of the i*^'^ dyon of 
kind m. Since the statistical weight of a config- 
uration of dyons is large when det G is large and 
small when it is small, det G imposes an attrac- 
tion between dyons that are nearest neighbours in 
kind, and a repulsion between same-kind dyons. 
The coefficients -1,2,-1 in front of the Coulomb 
interactions are actually the scalar products of 
the Cartan generators that determine the asymp- 
totics of the dyons' field, see subsection 6.2. 
The matrix G has the following nice properties: 



symmetry: G. 



rm,nj 



G. 



• overall "neutrality": the sum of Coulomb 
interactions in non-diagonal entries cancel 
those on the diagonal: Y,nj ^rni,nj = 4:TrUm 

• identity loss: dyons of the same kind are 
indistinguishable, meaning mathematically 
that det G is symmetric under permutation 
of any pair of dyons (i^j) of the same kind 
m. Dyons do not 'know' to which instanton 
they belong to 

• factorization: in the geometry when dyons 
fall into K well separated neutral clusters 
of N dyons of different kinds in each, det G 
factorizes into a product of exact integra- 
tion measures for K KvBLL instantons, 
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detG = (detG(i))^ where G^^^ is given by 
Eq. (63) 

• last but not least, the metric correspond- 
ing to G is hyper-Kahler, as it should be 
for the moduli space of a self-dual classical 
field [30]. In fact, it is a severe restriction 
on the metric. 

10. ENSEMBLE OF DYONS 

In the scmiclassical approximation we thus re- 
place the YM partition function (44) by the par- 
tition function of an interacting ensemble of an 
arbitrary number of dyons of N kinds: 

• detG(x„i), (68) 

where x^j is the coordinate of the i*^ dyon of 
kind TO, the matrix G is given by Eq. (67) and 
the fugacity / is given by Eq. (65). The overall 
exponent of the perturbative potential energy as 
function of the holonomy {/i} is understood, as 
in Eq. (64). 

The ensemble defined by a determinant of a 

matrix whose dimension is the number of par- 
ticles, is not a usual one. More customary, 
the interaction is given by the Boltzmann fac- 
tor exp (— ?7int(xi, . . .)). Of course, one can al- 
ways present the determinant in that way using 
the identity detG = exp(Tr logG) = exp(— f/mt) 
but the interactions will then include many-body 
forces. At the same time, it is precisely the de- 
terminant form of the interaction that makes the 
statistical physics of dyons an exactly solvable 
problem. 

10.1. Dyons' free energy: confining holon- 
omy preferred 

The partition function (68) can be computcxi 
directly and exactly, just by writing the deter- 
minant of G by definition as a sum of permuta- 
tions of products of the matrix entries. The result 
is astonishingly simple: all Coulomb interactions 
cancel exactly after integration over dyons' posi- 
tions, provided the overall neutrality condition is 



satisfied, viz. Ki = K2 = ■ ■ ■ = Kn = K; other- 
wise the partition function is divergent. We shall 
derive this interesting result by a more general 
method in section 11. Therefore, the recipe for 
computing the partition function is just to im- 
pose the neutrality condition and then to throw 
out all Coulomb interactions! We have thus to 
take the product of (47rj/m)'s from the diagonal 
of G: 

^ - 2^ ii '^^ ■ 

K ^ '' m=l 

The quantity AirfV is dimensionless and large for 
large volumes V. The sum can be therefore com- 
puted from the saddle point in K using the Stir- 
ling asymptotics for large factorials, and we ob- 
tain 

Z = exp {AnfVN {vii'2 . . . z/jv)"^) , 

1^1+1^2 + ■■■ + '^N = 1- (69) 

By definition, F = —TlogZ is the nonperturba- 

tivc dyon-induccd free energy as function of the 
holonomy; for = 2, 3 it is plotted in Fig. 9. 
Evidently, it has the minimum at 

1^1 = 1^2 = ■ ■ ■ = t^N = (70) 

corresponding to equidistant, that is confining 
values of fi's (47)! At the minimum, the free en- 
ergy is 

fmin = -Tl0gZ^in = -A7TfVT 

= -l^A^V^^^V (71) 

and there are no Coulomb corrections to this re- 
sult. In the last equation we have introduced the 
A'-indcpendcnt 't Hooft coupling A = asN/2'iT. 

We note that the free energy is 0{N'^) as ex- 
pected on general A''-counting grounds and that 
it is temperature-independent. It corresponds to 
\ogZ being proportional to the 4- volume V^'^^ = 
V/T, demonstrating the expected extensive be- 
haviour at low temperature. 

10.2. Cautionary remark on the numerical 
simulation of the dyon ensemble 
On the one hand, the determinant form for the 
weight of a dyon configuration is a very welcome 
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feature as it allows an exact computation not 
only of the free energy but also of various corre- 
lation functions, sec sections 11,12. On the other 
hand, it calls for questions that have been recently 
raised in Ref. [32]. The point is that the deter- 
minant (67) is not positive definite. Moreover, 
even if for some dyon configuration it is positive, 
individual eigenvalues of G need not be positive. 
If one takes many dyons on random in a dense 
ensemble, most of the eigenvalues will be nega- 
tive [32] . It rises the question of whether one can 
use the weight given by det G. 

Although the result for the free energy (69) 
looks as if the Coulomb interactions can be alto- 
gether neglected, the ensemble governed by det G 
is in fact very strongly correlated. The corre- 
lation functions computed in section 12 exhibit 
the Debye screening phenomenon typical for a 
Coulomb-like plasma. If one takes the eqiiidis- 
tant holonomy (70) minimizing the free energy 
one can present G — ^{E + C) where i? is a 
unity KN x KN matrix and C is a matrix formed 
by Coulomb interactions. One has then det G = 
const. exp(Tr log(l + C)) = const. exp(TrC — 
iTrC^ + . . .). As a matter of fact this weight is 
positive definite if the series is truncated at any 
finite number of terms. If one truncates it at the 
first term (which would be justifiable if the den- 
sity is low) it is precisely the Boltzmann weight 
for a Coulomb plasma. The Coulomb plasma, 
however, is not an exactly solvable system: there 
is an infinite series of corrections to the free en- 
ergy coming from interactions. The role of the 
higher terms in the expansion of log(l + C) is to 
cancel, order by order in density, those correc- 
tions to the free energy. The full series leads to 
a stronger correlation of dyons than even in the 
standard Coulomb plasma. This is clearly seen 
if one returns back to the det G weight: it goes 
to zero when two same-kind dyons come close 
enough together, corresponding to an infinitely 
high barrier in terms of the effective potential, 
whereas in the standard Coulomb plasma such 
configurations are only exponentially suppressed. 

Therefore, the ensemble governed by the det G 
weight strongly favours a clustering of dyons into 
neutral clusters of N different-kind dyons. For 
such configurations, all eigenvalues of G are posi- 



tive definite. In the limiting case of well-separated 
neutral clusters detG factorizes into a product 
of the one-instanton measures [27]. Therefore, if 
in a Monte Carlo simulation one starts with a 
well-clustered dyon configuration and then ther- 
molizes it by small enough Metropolis steps, the 
ensemble will never cross the line where a nega- 
tive eigenvalue appears, since, by continuity, the 
determinant must first become zero which would 
mean a zero weight for a configuration. However, 
implementing it in a practical simulation may be 
a diflacult task. 

11. STATISTICAL PHYSICS OF DYONS 
AS A QUANTUM FIELD THEORY 
IN 3D 

We follow here Ref. [33]. 

Although the Coulomb interactions of dyons 
cancel exactly in the free energy, the dyon en- 
semble defined by Eq. (68) is not a free gas but 
a highly correlated system. To facilitate comput- 
ing observables through correlation functions, we 
rewrite Eq. (68) as an equivalent quantum field 
theory. As a byproduct, we shall also check that 
the result for the free energy (69) is correct. 

To proceed to the quantum field theory de- 
scription we use two mathematical tricks. 

1. "Fermionization" (Berezin [34]). It is helpful 
to exponentiate the Coulomb interactions rather 
than keeping them in det G. To that end one 
presents the determinant of a matrix as an inte- 
gral over a finite number of anticommuting Grass- 
mann variables: 

det {Gab) = / H '^'^a d'lpA exp (^Va Gab ipB^ ■ 

Now we have the two-body Coulomb interac- 
tions in the exponent and it is possible to use 
the second trick presenting Coulomb interactions 
with the help of a functional integral over an 
auxiliary boson field. 

2. "Bosonization" (Polyakov [5]). One can 
write 
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= J D<p exp - Jdx (^'^^i^^i^ + P'pj 



An 

exp ( / p—p 



After applying the first trick the "charges" Q„i 
become Grassmann variables but after applying 
the second one, it becomes easy to integrate them 
out since the square of a Grassmann variable is 
zero. In fact one needs 2N boson fields Vm,Wm 
to reproduce diagonal elements of G and 2A'' an- 
ticommuting ("ghost") fields xL'Xm to present 
the non-diagonal elements. The chain of identi- 
ties is accomplished in Ref. [33] and the result for 
the partition function (68) is, identically, a path 
integral defining a quantum field theory in 3 di- 
mensions: 



Dx^ Dx Dv Dw exp / d'^x 



(72) 



— {dixlndiXm + diVmdiWm) 



+ f 

1 



{—^TTHm + Vm)^ hX^^ Xn 

dWm dWmOWn 



N 



(73) 



The subscript m is periodic: m = A'^ + 1 is equiva- 
lent to m = 1, and m = is equivalent to m = N. 

In the above path integral, the fields Vm have 
the meaning of the asymptotic Abelian electric 

potentials of dyons. 



(^4)„, 

A„4(x)/r = 27r/i„ - iwm(x). 



(74) 



while Wm have the meaning of the dual (or mag- 
netic) Abelian potentials. Note that the kinetic 
energy for the fields has only the mixing 

term diVmdiWm which is nothing but the Abelian 
duality transformation E • B. The function T{w) 
(73) is known as the periodic (or aflane) Toda lat- 
tice. 

Although the Lagrangian in Eq. (72) describes 
a highly nonlinear interacting quantum field the- 
ory, it is in fact exactly solvable! To prove it, one 
observes that the fields Vm enter the Lagrangian 
only linearly, therefore one can integrate them 



out. It leads to a functional (5-function: 
Dvjy,. — 



I 



This (5-function restricts possible fields Wm over 
which one still has to integrate in Eq. (72). 
Let Wm be a solution to the argument of the 
^-function. Integrating over small fluctuations 
about w gives the Jacobian 



Jac = det 



T 
47r' 



f 



dWmdWn 



(76) 



Remarkably, exactly the same functional determi- 
nant (but in the numerator) arises from integrat- 
ing over the ghost fields, for any background w. 
Therefore, all quantum corrections cancel exactly 
between the boson and ghost fields (a characteris- 
tic feature of supersymmetry), and the ensemble 
of dyons is basically governed by a classical field 
theory. 

This remarkable cancelation is of course not ac- 
cidental but can be traced to the hyper-Kahler 
nature of the dyons' metric, which in its turn 
stems from the self-duality of dyons. 

To find the ground state we examine the fields' 
potential energy being —AnfiJmd^/dwm, which 
we prefer to write restoring Vm = Mm+i — Mm and 

as 



(77) 



(the volume factor arises for constant fields Wm)- 
One has first to find the stationary point in Wm 
for a given set of fm's- It leads to the equations 



dV 

dWm. 







whose solution is 

„w,-W2 (l/ll/2l^3-i'iv)^ 



etc. 



Putting it back into Eq. (77) we obtain 

r = -ATTfVN{viV'2--l'N)^ , 
Vl+l'2 + ■■■ + !^JV = 1, 



(78) 
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which is exactly what one gets from a direct cal- 
culation of the partition function, outlined in the 
previous section, sec Eq. (69). The minimum is 
achieved at the equidistant, confining value of the 
holonomy, see Eqs.(70, 47). Using field-theoretic 
methods, we have also proven that the result is 
exact, as all potential quantum corrections can- 
cel. It is in line with the exact cancelation of the 
Coulomb interactions in the determinant. 

Given this cancelation, the key finding - that 
the dyon-induced free energy has the minimum 
at the confining value of holonomy - is in a sense 
trivial. If all Coulomb interactions cancel af- 
ter integration over dyons' positions, the weight 
of a many-dyon configuration is the same as if 
they were infinitely dilute (although they are). 
Then the weight, what concerns the holonomy, 
is proportional to the product of diagonal ma- 
trix elements of G in the dilute limit, that is to 
the normalization integrals for dyon zero modes. 
These arc nothing but the field strengths i*],,^ 
of individual dyons, hence the normalization is 
proportional to the product of the dyon actions 
Sm = STTt'm/as. The sum of all TV kinds of dyons' 
actions is fixed and equal to the instanton action, 
however, it is the product of actions that defines 
the weight. The product is maximal wlic^n all ac- 
tions are equal, hence the equidistant or confining 
/u's are statistically preferred. Thus, the average 
Polyakov line is zero, < Tri >= 0. 

12. HEAVY QUARK POTENTIAL 

The field-theoretic representation of the dyon 
ensemble enables one to compute various YM cor- 
relation functions in the scmiclassical approxima- 
tion. The key observables relevant to confinement 
are the correlation function of two Polyakov lines 
(defining the heavy quark potential), and the av- 
erage of large Wilson loops. A detailed calcula- 
tion of these quantities is performed in Ref . [33] ; 
here we only present the results and discuss the 
meaning. 

12.1. A^-ality and fc-strings 

From the viewpoint of confinement, all irre- 
ducible representations of the SU{N) group fall 
into N classes: those that appear in the direct 



product of any number of adjoint representations, 
and those that appear in the direct product of 
any number of adjoint representations with the 
irreducible representation being the rank-A; anti- 
symmetric tensor, k = 1, . . . ,N —1. 'W-ality" 
is said to be zero in the first case and equal to k 
in the second. A''-ality-zero representations trans- 
form trivially under the center of the group Z^; 
the rest acquire a phase 2'Kk/N . 

One expects that there is no asymptotic linear 
potential between static colour sources in the ad- 
joint representation as such sources are screened 
by gluons. If a representation is found in a di- 
rect product of some number of adjoint represen- 
tations and a rank-fc antisymmetric representa- 
tion, the adjoint ones "cancel out" as they can 
be all screened by an appropriate number of glu- 
ons. Therefore, from the confinement viewpoint 
all A''-ality = k representations are equivalent and 
there are only N — 1 string tensions ak,N being 
the coefficients in the asymptotic linear potential 
for sources in the antisymmetric rank-fc represen- 
tation. They are called "fc-strings". The repre- 
sentation dimension is dk,N = k\{N-ky. ^^'^ 
eigenvalue of the quadratic Casimir operator is 
Ck,N = ^HN-k). 

The value k = l corresponds to the fundamental 
representation whereas k = N—1 corresponds to 
the representation conjugate to the fundamental 
[quarks and anti-quarks]. In general, the rank- 
{N—k) antisymmetric representation is conjugate 
to the rank-fc one; it has the same dimension and 
the same string tension, ak,N = o^N-k,N■ There- 
fore, for odd all string tensions appear in equal 
pairs; for even N, apart from pairs, there is one 
privileged representation with fc = ^ which has 
no pair and is real. The total number of different 
string tensions is thus [^] . 

The behaviour of ak,N as function of fc and N 
is an important issue as it discriminates between 
various confinement mechanisms. On general N- 
counting grounds one can only infer that at large 
AT and fc < N, ak,N/<Ji,N = {k/N){l+0{l/N^)). 
Important, there should be no 0{N~^) correc- 
tion [35] . A popular version called "Casimir scal- 
ing" , according to which the string tension is pro- 
portional to the Casimir operator for a given rep- 
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resentation (it stems from an idea that confine- 
ment is somehow related to the modification of a 
one-gluon exchange at large distances), does not 
satisfy this restriction 

12.2. Correlation function of Polycikov 

lines 

To find the potential energy Vk,N of static 
"quark" and "antiquark" transforming according 
to the antisymmetric rank-fc representation, one 
has to consider the correlation of Polyakov lines 
in the appropriate representation: 



TrLk,N{zi) TV4,Jv(z2) 

Vfc,jv(zi - Z2) 



(79) 



const, exp 



Far away from dyons' cores the field is Abelian 
and in the field-theoretic language of Eq. (72) is 

given by Eq. (74). Therefore, the Polyakov line 
in the fundamental representation is 



N 



TVL(z) = J2^" 



(80) 



m=l 



exp I lirifjin 



2^™ (=2) 



In the general antisymmetric rank-A; representa- 
tion 

N 

TrLfe,jv(z)= X! 

(81) 

mi<m2<...<mfc 

where cyclic summation from 1 to iV is assumed. 

The average (79) can be computed from the 
quantum field theory (72). Inserting the two 
Polyakov lines (81) into Eq. (72) we observe that 
the Abelian electric potential Vm enters linearly 
in the exponent as before. Therefore, it can be in- 
tegrated out, leading to a 5-function for the dual 
field Wm, which is now shifted by the source (cf. 
Eq. (75)): 



Dv„ 



T 
Ait 



+ / 



dWrr,. 



^There are many lattice studies claiming Casimir scaling 
for various groups and representations. It refers to certain 

"intermediate" separations between static sources. We are 
now discussing the string tension at asymptotic separa- 
tions. 



- -(5(x-Zi)(5„ 



mi + • • • + Smmk 



(82) 



One has to find the dual field u'm(x) nullifying 
the argument of this 5-function, plug it into the 

action 



exp dx ^^(w)^ > 



(83) 



and sum over all sets {mi < m2 < ... < ruk}, 
{rii < n2 < ... < rifc} with the weight 
exp (27ri(mi -I- . . . -f nik — ni — ... — nk)/N). 
The Jacobian from resolving the (5-function again 
cancels exactly with the determinant arising from 
ghosts. Therefore, the calculation of the correla- 
tor (79), sketched above, is exact. 

At large separations between the sources |zi — 
Z2I, the fields w„i resolving the (5-function are 
small and one can expand the Toda chain: 



dWrn 

where 
M = 



N+-Wr, 



Mr 



,Wr, 



Mr 



(84) 
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(85) 



\ -1 ... -1 2 / 

As apparent from Eq. (84) , the eigenvalues of M 
determine the spectrum of the dual fields Wm- 
There is one zero eigenvalue which decouples from 
everywhere, and N—1 nonzero eigenvalues 



nk 
2sm- 



1,...,N- 



1. 



(86) 



Certain orthogonality relation imposes the selec- 
tion rule: the asymptotics of the correlation func- 
tion of two Polyakov lines in the antisymmetric 
rank-fc representation is determined by precisely 
the k*"^ eigenvalue. We obtain [33] 



TrLfc,jv(zi) Tr4_j^(z2)^ 

°° const, exp (-|zi - zaj M^A^ ^(87) 
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where M is the 'dual photon' mass, 



'neighbours in kind' are connected [21]: 



M = 



~Y~ ~ 2tt\T 



0{N). 



(88) 



Comparing Eq. (87) with the definition of the 
heavy quark potential (79) we find that there is 

an asymptotically linear potential between static 
"quarks" in any iV-ality nonzero representation, 
with the fc-string tension 



(jfc^jv = MtVxW =2MT sin^ 

A2 iV . Trfc 
= sm — 

A TT N 



(89) 



independent of temperature and of N at large A'' 
and fixed k. This is the so-called 'sine regime': it 
has been found before in certain supersymmetric 
theories [36] . Lattice simulations [37] support this 
regime, whereas another lattice study [38] gives 
somewhat smaller values but within two stan- 
dard deviations from the values following from 
Eq. (89). For a general discussion of the sine 
regime for A:-strings, which is favoured from many 
viewpoints, see [35]. 

We see that at large N and k <C N , 
crk,N/cri,N = {k/N){l + 0{l/N^)), as it should 
be on general grounds, and that all fc-string ten- 
sions have a finite limit at zero temperature. 

12.3. Debye screening and clustering of 
dyons 

Essentially the same calculation gives the corre- 
lation function of dyons in the ensemble defined 
by Eq. (68). The correlation is large at small 
distances and decays exponentially at large dis- 
tances, see Eq. (87). Such behaviour is known as 
Debye screening and is typical for the Coulomb 
plasma. Screening implies that every dyon of any 
kind is on the average surrounded by iV — 1 dyons 
of other kinds at a range of the order or less than 
the Debye radius ro ~ 1/M. Therefore, the dyon 
ensemble is in fact strongly correlated into neu- 
tral clusters, as discussed in subsection 10.2. 

A neutral cluster of N different kinds of dyons 
is in fact an instanton (section 7). It is interesting 
to estimate its average size. The KvBLL instan- 
ton size is related to the perimeter of the polygon 
formed by N constituent dyons where the nearest 
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' m,m+l 

27rT ' 



(90) 



Assuming for an estimate that the average separa- 
tion between dyons of neighbour kind in a cluster 
is of the order of the Debye radius, rm,m+i ^ fo, 
we find that the average instanton size is 



<p>' 



A 



(91) 



being independent either of temperature or oi N . 
This is the expected typical size of the Ad lumps 
formed by correlated neutral clusters of N dyons. 

Let us compare the Debye radius Td ''^ ^ = 
^J^ , and the size of the core of an individual 
dyon Tcore ~ 2-kTv ~ ' subscctiou 6.2. We 
observe that at temperatures T < the De- 

bye range lies within the cores of individual dyons. 
Since the deconfinement temperature is expected 
to be Tc A which is less than the above esti- 
mate, it means that in the confinement phase the 
screening occurs inside the typical core range at 
all temperatures! This may seem paradoxical but 
one has to remember that the measure (63) as due 
to zero modes is exact even if dyons overlap, and 
that it is this measure that leads to the Debye- 
like screening and to the clustering of dyons into 
neutral clusters. 

13. AREA LAW FOR LARGE WILSON 
LOOPS 

When dealing with the ensemble of dyons, it 
is convenient to use a gauge where Ai is diago- 
nal {i.e. Abelian). This necessarily implies Dirac 
string singularities sticking from dyons, which are 
however gauge artifacts as they do not carry any 
energy, see subsection 6.1. Moreover, the Dirac 
strings' directions are also subject to the freedom 
of the gauge choice. For example, one can choose 
the gauge in which N dyons belonging to a neu- 
tral cluster are connected by Dirac strings. This 
choice is, however, not convenient for the ensem- 
ble as dyons have to loose their "memory" to what 
particular instanton they belong to. The natural 
gauge is where all Dirac strings of all dyons are 
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directed to infinity along some axis, e.g. along 
the z axis, see subsection 6.1. 

In this gauge, the magnetic field of dyons be- 
yond their cores is Abelian and is a superposition 
of the Abelian fields of individual dyons. For large 
Wilson loops we are interested in, it is this su- 
perposition field of a large number of dyons that 
contributes most as they have a slowly decreas- 
ing l/|x — Xi| asymptotics, hence the use of the 
field outside the cores is justified. Owing to self- 
duality, 

T 

-- --^ i5m„5iWm(x),(92) 



[i?,(x)],„„ = [d,A,{^)i 



cf. Eq. (74). Since Ai is Abelian beyond the 
cores, one can use the Stokes theorem for the spa- 
tial Wilson loop: 

W = TiVexpi j>Aidx' = Tr exp i j Bi dV' 

= (-4 / "^'^^ ^'""A • ^^^^ 

Eq. (93) may look contradictory as we first use 
Bi = cur Mi and then Bi = diA^. Actually there 
is no contradiction as the last equation is true up 
to Dirac string singularities which carry away the 
magnetic flux. If the Dirac string pierces the sur- 
face spanning the loop it gives a quantized con- 
tribution exp(27rz- integer) = 1; one can also use 
the gauge freedom to direct Dirac strings paral- 
lel to the loop surface in which case there is no 
contribution from the Dirac strings at all. 

Let us take a flat Wilson loop lying in the [xy) 
plane at z = 0. Then Eq. (93) is continued as 



W 



2^ oxp I -I 
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(fx d:^VmS{z) 
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exp I I— 



a:,yGArea 
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It means that the average of the Wilson loop in 
the dyons ensemble is given by the partition func- 
tion (72) with the source 

E / ^™ y ^ Area) j 

where 0{x,y G Area) is a step function equal to 
unity if x,y belong to the area inside the loop 



and equal to zero otherwise. As in the case of the 
Polyakov lines the presence of the Wilson loop 
shifts the argument of the ^-function arising from 
the integration over the Vm variables, and the av- 
erage Wilson loop in the fundamental representa- 
tion is given by the equation 



{W) 
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mi 



47r/ 

Dwm exp ( / dx — -^(w^) 



■ det --j-a Smn + ^ — 



n 



iT 
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dS{z) 



- 9{x, y G Area) J .(95) 



Therefore, one has to solve the non-linear equa- 
tions on Wm's with a source on the surface of the 
loop, 

dS{z 



-2TTi (5. 



mmi —J — dix, y e Area), (96) 
dz 



for all mi, plug it into the action {^-k f / N)T (w) , 
and sum over m\. In order to evaluate the aver- 
age of the Wilson loop in a general antisymmetric 
rank-A: representation, one has to take the source 
in Eq. (96) as -2Tii5'{z) ((5mmi+- • ■+^mmk) and 
sum over mi < . . . < rrik from 1 to A'', see Eq. (81). 
Again, the ghost determinant cancels exactly the 
Jacobian from the fluctuations of Wm about the 
solution, therefore the classical-field calculation is 
exact. 

Contrary to the case of the Polyakov lines, one 
cannot, generally speaking, linearize Eq. (96) in 
Wm but has to solve the non-linear equations as 
they are. The Toda equations (96) with a 5'{z) 
source in the r.h.s. define "pinned soliton' solu- 
tions Wm{z) that are Id functions in the direc- 
tion transverse to the surface spanning the Wil- 
son loop but do not depend on the coordinates 
X, y provided they are taken inside the loop. Be- 
yond that surface Wm = 0. Along the perimeter of 
the loop, Wm interpolate between the soliton and 
zero. For large areas, the action (83) is therefore 
proportional to the area of the surface spanning 
the loop, which gives the famous area law for the 
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average Wilson loop. The coefBcient in the area 
law, the 'magnetic' string tension, is found from 
integrating the action density of the solution in 
the z direction. 

The exact solutions of Eq. (96) for any N and 
any representation k have been found in Rcf . [33] , 
and the resulting 'magnetic' string tension turns 
out to be 



N 



Trk 



(^k,N = sm — , 

An TV 



(97) 



which coincides with the 'electric' string tension 
(89) found from the correlators of the Polyakov 
lines, for all A;-strings! 

Several comments are in order here. 

• The 'electric' and 'magnetic' string tensions 

should coincide only in the limit T — !■ 
where the Euclidean 0(4) symmetry is re- 
stored. Both calculations have been in 
fact performed in that limit as we have 
ignored the temperature-dependent pertur- 
bative potential (49). If it is included, 
the 'electric' and 'magnetic' string tensions 
split. 

• despite that the theory (72) is 3- 
dimensional, with the temperature entering 
just as a parameter in the Lagrangian, it 
"knows" about the restoration of Euclidean 
0(4) symmetry at T — » 0. 

• the 'electric' and 'magnetic' string tensions 
are technically obtained in very different 
ways: the first is related to the mass of the 
elementary excitation of the dual fields Wm, 
whereas the latter is related to the mass of 
the dual field soliton. 

14. ADDING ANTI-DYONS 

So far we have considered the YM vacuum filled 
with self-dual dyon solutions only. However, the 
CP invariance of the vacuum requests, in the 
language of the semiclassical approximation, that 
there is an equal number of anti-self-dual dyons in 
the vacuum, up to thermodynamic fluctuations of 
the order of VK \/V. A superposition of self- 
dual and anti-self-dual solutions is not, strictly 



speaking, a solution of the YM equation of mo- 
tion anymore, unless they are well separated. 

In particular, at distances exceeding the core 
sizes, dyons and anti-dyons experience classical 
Coulomb interactions (see e.g. Ref. [39]) 
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where x^j is the position of the i^^ dyon of kind 
m, and ymi is the position of the i*'^ anti-dyon of 
kind m. We expect, however, that at separations 
larger than the inverse Debye screening mass 1/M 
which in fact is smaller than the core size (see 
subsection 12.3) the classical Coulomb interaction 
is totally screened and therefore can be neglected. 

Dyons and anti-dyons experience also quantum 
interactions that are of the order of unity in the 
coupling constant g^. These arise from («) non- 
factorization of the fluctuation determinant, {ii) 
non-factorization of the integration measure. It is 
quite difficult to define these quantities properly 
as the dyon-anti-dyon configuration is not a strict 
saddle point. The topic has been much discussed 
in the context of building the 'instanton liquid' 
model of the YM vacuum where a similar prob- 
lem arises with respect to instantons and anti- 
instantons [13]. There is a way to define properly 
both 'classical' and 'quantum' interactions of in- 
stantons and anti-instantons by means of unitar- 
ity and analyticity, however for dyons this ap- 
proach has not been developed. 

Therefore, we adopt here a simple model that 
the leading source of the dyon anti-dyon interac- 
tion is (ii). By definition, the integration measure 
over the collective coordinates is the determinant 
of a matrix built of the overlaps of the (quasi) 
zero modes. It is known exactly for a neutral 
cluster of N kinds of dyons, the KvBLL caloron, 
see Eq. (63), and approximately for many dyons, 
Eq. (67). The {K+N x K+N) matrix G, where 
K+ is the (equal) number of dyons of any kind m, 
will now become a {{K+ + K^)Nx {K+ + K^)N) 
matrix Q where the diagonal {K±N x K±N) 
blocks are built from the overlaps for the same- 
duality dyons whereas the off-diagonal {K±N x 
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K^N) blocks describe the interaction of opposite- 
duality dyons, that is dyons and K_ anti- 
dyons. 

The diagonal blocks for same-duality dyons are 
the same as before (sec (67)), with the only mod- 
ification that the diagonal elements composed of 
the normalization integrals for zero modes, being 
the dyon action and thus proportional to A4, get 
Coulomb corrections also from opposite-duality 
dyons. Therefore, the first two lines in Eq. (67) 
acquire an addition 



mk I 



^p l_ 



^ T\'X.mi — ym+l,k\ ^ T\^mi — ym-\.k\ j 

The off-diagonal blocks for opposite-duality 
dyons are composed of Coulomb bonds as in 
the second two lines of Eq. (67), with the re- 



placement |x„ 



Note 



that the signs of all Coulomb bonds, both at 
the diagonal and off-diagonal, are the same as in 
the case of same-duality overlaps. The resulting 
+ KJ)N X {K+ -h K_)N) matrix Q whose 
determinant defines the dyon-anti-dyon ensemble; 
retains the important property that all Coulomb 
bonds cancel in the sum of the elements of Q along 
any row or any column. 

With this natural modification of the measure 
over the moduli space to incorporate anti-dyons, 
one can again rewrite the partition function of 
the ensemble, now composed of dyons and anti- 
dyons, as a path integral for a 2>d quantum field 
theory. It is exactly the same theory as given by 
Eq. (72) but with a double fugacity, / 2/. The 
derivation of the 'electric' and 'magnetic' string 
tensions goes as before (see sections 12 and 13) 
with the simple replacement a \/2a. 

We note that had we altogether neglected the 
interactions between dyons and anti-dyons treat- 
ing them as two independent 'liquids' we would 
have to introduce two sets of auxiliary fields 
w, w, x^, X to write down the ensemble as a quan- 
tum field theory. The free energy of the system 
would then also double but the correlation func- 
tions would be different. Meanwhile, the boson 



fields Vm have the meaning of the Abelian po- 
tentials (74) and w^'s have the meaning of the 
dual Abelian potentials, and there should be no 
'double sets' of those fields in the effective action. 
This is achieved naturally if one introduces the 
dyon-anti-dyon interaction through the overlaps 
of the quasi-zero modes as described above. 

The proposed model for accounting anti-dyons 
is still a guesswork inviting for a better under- 
standing. Nevertheless, a quantity that is sensi- 
tive to the dyon-anti-dyon interactions - the so- 
called topological susceptibility - turns out to be 
quite reasonable. 

14.1. Topological susceptibility and the 
gluon condensate 

An important characterization of the vacuum 
is the topological susceptibility defined as the 
zero-momentum correlation function of topologi- 
cal charge densities: 



hiFF{x) TrFF(0)\ 



167r2 



167r2 



If one adds the 6 term to the action. 



i0 
167r2 



/ 



rJ'^rTr F F 



and computes the free energy F of the system, 
the topological susceptibility can be obtained as 



?T> = 



1 d^F 



V 9612 



61=0 



The m dyon carries the topological charge i^„j 
whereas the m**^ anti-dyon carries the topological 
charge —Vm, see Eq. (58). Therefore from the 
point of view of dyons, the introduction of the 6 
term changes the fugacity of the w}^ dyon / — > 
/ e'^'^"' and the fugacity of the m*^ anti-dyon / — > 
/ e~*^'^" . The potential energy of the ensemble 
(77) changes therefore to 



Minimizing it in and Vm we find again that 
the confining holonomy i/„ = -i- gives the mini- 
mum, and the free energy at the minimum is (cf. 
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Eq. (71)) 



F = -Stt/ VT cos — = 



Differentiating it twice with respect to 6 we find 
the topological susceptibility 

27r2 A2 



<Qt>= 7T-^TT 



We see that the topological susceptibility is stable 
at large N as it is expected from the A''-counting 
rules. 

Another important quantity characterizing the 
vacuum is the gluon condensate also related, via 
the trace anomaly, to the free energy (see e.g. 
Ref. [13]): 



<TrJ-^^.> 
167r2 
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= N- 
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A4 



> 
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The A''-dependence of the condensate is also as 

expected. 

Let us calculate the average density of dyons 
(and anti-dyons) in the ensemble. To that end, 
we artificially introduce separate fugacities for 
dyons of the m^^ kind into the partition function 
(68); then the average number of dyons is found 
from the obvious relation 



<Km> = 



dlogZ 



dlogfn 



fm=f 



With separate fugacities, the result (71) is modi- 
fied by replacing / (/i/2---/jv)^/^, hence 

dF 



<Kr, 



^ rp fn 



dfrr 
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i.e. a finite and equal density of each kind of 
dyons and anti-dyons in the 4- volume F , mean- 
ing also the finite density of the KvBLL instan- 
tons and anti-instantons. From the 3-dimensional 
point of view, the 3d density of dyons (and 
KvBLL instantons) is increasing as the tempera- 
ture goes down: there are more and more instan- 
tons sitting on top of each other in 3rf but spread 
over the (compactified) time direction. 



Although we have not expected that our theory 
is valid at small temperatures (where the approxi- 
mate measure we use for same- kind and opposite- 
duality dyons is probably incomplete since dyons 
are dense) the zero-temperature limit is obtained 
in a smooth way. We have seen it already exam- 
ining the 'electric' and 'magnetic' string tensions, 
and now we see that the topological susceptibility 
and the gluon condensate have also a reasonable 
zero-temperature limit, while the logarithm of the 
partition function is extensive in the Ad volume 
as it should be. It is interesting that these re- 
sults are found from a theory that is essentially 
3-dimensional. A popular word used in such sit- 
uation is "holography" : the 3c? theory apparently 
'knows' about and adequately describes the 4d 
world. 

Finally, let us make a numerical estimate of the 
calculated quantities. For a numerical estimate 
at N = 3, we take A = 1/4 compatible with the 
commonly assumed freezing of as at the value 
of 0.5, and A = 200 MeV in the Pauli-Villars 
scheme. We then obtain the topological suscep- 
tibility, the gluon condensate and the string ten- 
sion (190MeV)4, (255MeV)4 and (433 MeV)^, 
respectively, being in good agreement with the 
phenomenological and lattice values. 



15. CANCELATION OF GLUONS IN 
THE CONFINEMENT PHASE 

To prove confinement, it is insufficient to 
demonstrate the area law for large Wilson loops 
and the zero average for the Polyakov line: it 
must be shown that there are no massless glu- 
ons left in the spectrum. We give an argument 
that this indeed happens in the dyon vacuum. 

A manifestation of massless gluons in perturba- 
tion theory is the Stefan-Boltzmann law for the 
free energy: 



FsB T 

— = --l0g^SB 



: 2 (number of gluons) J ^ log (l ^ e j 
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(99) 



It is proportional to the number of gluons N'^ — 1 
and has the T"* behaviour characteristic of mass- 
less particles. In the confinement phase, neither 
is permissible: If only glueballs are left in the 
spectrum the free energy must be 0{N^) and the 
temperature dependence must be very weak un- 
til T Pi Tc where it abruptly rises owing to the 
excitation of many glueballs. 

As explained in Section 8, the ensemble of 
dyons has a nonperturbative free energy 



7V2 
■27r2 A2" 



(100) 



-^dyon 

V ~ 

It is 0{N^) but temperature-independent. Dyons 
force the system to have the 'most nontrivial' 
holonomy (47). For that holonomy, the pertur- 
bative potential energy (49) is at its maximum 
equal to 



■ pert, max 

V 



45 
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(101) 



The full free energy is the sum of the three terms 
above: 



1 

7V2 



(102) 
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We see that the leading 0{N'^) term in the 
Stefan-Boltzmann law is canceled by the poten- 
tial energy. It happens only at precisely the 'most 
nontrivial', confining value of the holonomy and 
nowhere else! In fact it seems to be the only way 
how 0{N'^) massless gluons can be canceled out 
of the free energy, and the main question shifts to 
why docs the system prefer the 'confining' holon- 
omy. Dyons seem to be in a position to answer 
the question. 

There is a 0{T^N^) term left in Eq. (102), 
which is also unacceptable in the confinement 
phase and has to be canceled, too. However at 
the 0{N'^) level there appear also quantmn cor- 
rections to the free energy not considered here. 
Although it has not been checked explicitly there 
is a good chance that all O(T^) terms in the free 



energy will be absent since there are no massless 
degrees of freedom left in the theory. Actually 
the spectrum is determined by string excitations 
but they are not studied yet. 

16. DECONFINEMENT PHASE TRAN- 
SITION 

As the temperature rises, the perturbative free 
energy grows as and eventually it overcomes 
the negative nonperturbative free energy (100). 
At this point, the trivial holonomy for which both 
the perturbative and nonperturbative free energy 
are zero, becomes favourable. Therefore an es- 
timate of the critical deconfinement temperature 
comes from equating the sum of Eq. (100) and 
Eq. (101) to zero, which gives 

(103) 



^5 
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As expected, it is stable in A^. A more robust 
quantity, both from the theoretical and lattice 

viewpoints, is the ratio T^/ ^/a since in this ra- 
tio the poorly known parameters A and A cancel 
out: 
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(104) 





su{i) 


SU{A) 


SU{Q) 


SU{S) 


theory 


0.6430 


0.6150 


0.5967 


0.5906 


lattice 


0.6462 
(30) 


0.6344 
(81) 


0.6101 
(51) 


0.5928 
(107) 



In the Table, we compare the values from 
Eq. (104) to those measured in lattice simulations 
of the pure SU {N) gauge theories [40] ; there is a 
surprisingly good agreement. 
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For the SU{2) gauge group, the phase transi- 
tion is of the 2nd order, as seen from Fig. 14; for 
higher-rank groups it is of the 1st order. 



potential energy in SU{2) 




Figure 14. Potential energy of the pure SU{2) 
YM theory as function of ^TrL at zero {lower 
curve), critical [middle curve) and high [upper 
curve) temperatures. At the critical temperature 
the curve is flat exhibiting the second order phase 
transition. 



To illustrate the confincment-deconfinement 
phase transition, we sum the nonperturbative po- 
tential energy following from (69) and the pertur- 
bative potential energy (49), and plot it for the 
5t/(3) gauge group as function of the holonomy 
in Fig. 15. The contour plot has the symmetry of 
the extended root diagram for the group and is 
double periodic in A\ and A\, if one takes the di- 
agonal components of as a basis. One can con- 
centrate on the fundamental domain correspond- 
ing to the Weyl camera in the root diagram; all 
the rest are obtained from this one by shifts and 
reflections. At zero temperature the perturbative 
potential energy is zero while the nonperturba- 
tive one plotted in the lower part of Fig. 9 has 
the minimum at the 'confining holonomy' corre- 
sponding to TrL = 0. It is the confining phase. 

As the temperature rises, the perturbative po- 
tential energy starts to play an increasingly im- 
portant role. At high temperature it is dominant, 
and asymptotically the minimum is at the trivial 
field A\ = A\ = Q (at the center of the plots), 




Figure 15. Contour plots for the potential energy 
in the (Al,^!!) plane at zero [upper left), critical 
[upper right) and high [lower left) temperatures 
for the SU[i) gauge group. A zoom at critical 
temperature [lower right) clearly shows the coex- 
istence of two minima. 



and at its copies by periodicity. At a temper- 
ature close to that estimated in (103) there are 
two minima: one is at the 'confining' holonomy, 
the other is at small values of A\, A\. Such a 
situation is referred to as the phase transition of 
the 1st order: when one of the two minima be- 
comes more deep than the other the system jumps 
to the deeper one; the other is then meta-stable. 
The first derivative of the free energy e.g. with 
respect to the temperature is discontinuous as the 
system jumps. 

17. CONFINEMENT AND DECON- 
FINEMENT IN THE G(2) GROUP 

It is interesting to check the dyon mechanism 
of confinement and of deconfinement in the G(2) 
gauge group. G(2) is an exceptional group of 
rank 2 with 14 generators, hence 14 gluons, out of 
which 8 are the same as in SU[2>) and the other 
6 have the colour quantum numbers of 3 quarks 
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and 3 anti-quarks. 

G(2) can be considered as a subgroup of 5*0(7) 
whose group elements are represented by 7 x 7 
orthogonal matrices (QQ^ = 1, detfi = 1) with 
an additional constraint 

where T^bc is an antisymmetric matrix with the 
following nonzero elements: 

1^127 = Ti54 = Ti63 = T235 

= 7264 = 7374 = T576 = 1, TacdTbcd = Q^ab- 

The lowest-dimensional representation of G{2) is 
7 (we shall call it fundamental), the next are 14 
(the adjoint representation), 27, etc. 

The main reason why G(2) gauge theory is in- 
teresting is that the G{2) group has a trivial cen- 
ter: only the unity matrix commutes with all the 
other group elements. In contrast, in the SU{N) 
group the center is nontrivial since N group ele- 
ments of the form exp(27rifc/iV)ljv, k = 1,...N, 
commute with other elements of the group; they 
form a discrete subgroup of SU{N). There- 
fore, the G(2) gauge theory serves as an excel- 
lent laboratory to check various ideas about the 
mcichanism of confinement. In particular, it is 
almost universally accepted that confinement in 
the SU{N) gauge theory is related to the cen- 
ter group Zn, and that deconfinement is directly 
related to the spontaneous breaking of the sym- 
metry with respect to Zjv. If the mechanism of 
confinement is universal for all Lie groups and if 
there is confinement in the centerless G(2) group, 
this popular view has to be reconsidered. 

The gauge theory based on the G(2) group has 
been extensively studied on the lattice in the last 
few years by several groups [41]- [44]. The conclu- 
sion is unanimous: There is confinement at low 
temperature, in the sense that < Tr L >= where 
L is the Polyakov line in the fundamental repre- 
sentation 7, and there is a Ist-order deconfine- 
ment transition to a phase where <TvL>^ 0, - 
just as in the SU (3) gauge theory. It means that 
the presence or the absence of a nontrivial center 
of the gauge group is irrelevant to the mecha- 
nism of confinement, and that the deconfinement 
transition is unrelated to the breaking of center 



symmetry, as there is no such symmetry in G(2). 
As a matter of fact this has been long advocated 
by Smilga [45] on rather general grounds. 

The question is whether dyons are capable to 
explain confinement and deconfinement in the 
G(2) gauge theory, as they manage to do it for 
SU{N). The answer is "yes", and it serves as 
a very nontrivial check of the philosophy and of 
the formalism. The material below is based on 
the work in collaboration with Victor Petrov, in 
preparation. 

17.1. Dyons in the G(2) gauge group 

There are 14 generators in G(2); in the lowest 
7-dimensional representation eight of them can be 
chosen in the form 



A„ 




Tr (AaAa) = 1, 



where Aq are the standard SU{3) Gell-Mann ma- 
trices. We choose the two Cartan generators 
H^ g, = A^^g in similarity with SU{3). There are 
12 "shift" generators Ef satisfying the commu- 
tation relations 

[H3, Ef] = a^iEf, [Hs, Ef] = a^^Ef . 

The set of 12 two-dimensional vectors cXi = 
{oi^i,cts,i) forms a root system of G2. We 
choose ai = (—5)""^)) cti = 1) and 

012 = (5)^^) I c^i = 51 a^s the simple roots. 

All the rest roots are linear combinations of 
these. We also introduce the third root cxq= 
^— ^, , olq — 1, which is not simple but has 

the largest negative coefficients when expanded 
in simple roots. Correspondingly, there are three 
types of fundamental dyons in G(2), whose mag- 
netic charges are determined by the dual roots. 



mo = - 



mi 



2' 2 J 

1 _vi' 

2' 2 



(105) 
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It implies that the asymptotic electric and mag- 
netic fields of those three dyons are 



E, 



0,1,2 



Bo,i,2 = -j— To (mo, 1,2 • H) 



(106) 



where H = (Hs.Hfi) is a diagonal matrix. The 
dual roots appear here because the dyon is an 
SU{2) object: When constructing it along the 
lines of section 6 one picks up an SU (2) subgroup 
whose generators commuting as the three Pauli 
matrices are 



72. 



E+-E7 



ICXi 



As usually, dyon solutions are characterized by 
the holonomy or the eigenvalues of the Polyakov 
line (here: in the 7-dimcnsional representation) at 
spatial infinity. In the gauge where A4 is static 
and diagonal we label the holonomy by two num- 
bers fi = (hsjIIs) such that asymptotically 

^4(00) = 2^T(/x-H) = 27:T{tX3H3 + nsHs).{l07) 

Correspondingly, the Polyakov line's eigenvalues 
are 



L = cxp (27ri(/i, • H)) 



(108) 



The actions of the three individual dyons are 

5i,2 = —(/X -1111,2), (109) 
as 

Ait 

So = — (l + /x-mo). 

as 

The above construction (106)-(110) can be gen- 
eralized to an arbitrary gauge group [17]: One 
takes r dyon solutions associated with the simple 
roots ai, i = l,...,r, where r is the rank, and 
supplements them by a dyon associated with a 
non-simple lowest root cxo that has the most neg- 
ative coefficients in the expansion over the simple 
roots. In the gauge where A4 is time indepen- 
dent the first r dyons are static while the last 
or zeroth dyon is obtained by a time-dependent 
gauge transformation; in the notations of subsec- 
tion 6.1 it is the "L" dyon. The magnetic charges 
of all r -|- 1 dyons are given by the dual roots: 



m,; 



0,l,...,r. 



Since mo is not linearly independent as an r- 
dimensional vector it can be expanded in the 
other dual roots, mo = — K^i^ 



i=0 



k*mi = 0, 



(110) 



where the integers k* are called dual Kac labels, 
or co-marks; by construction fco = 1. In the case 
of G{2) we have fcg = 1, kl = 2, k^ = 1. The 
sum 



C2 



i=0 



(111) 



(equal to 4 for G(2) and N for SU{N)) is called 
the dual Coxeter number. Eq. (110) tells us that 
in order to build an electric- and magnetic-neutral 
object (the KvBLL instanton with a nontriv- 
ial holonomy) one needs C2 fundamental dyons, 
some of which may have a multiplicity other than 
unity. For example, in the G{2) gauge theory one 
has to take the dyon of the 1st kind twice since 
mo -|- 2mi -|- m2 = 0, meaning that in this case 
the KvBLL instanton is made of four dyons but 
only of three kinds. As a curiosity, the KvBLL 
instanton of the exceptional E{8) gauge theory is 
built out of 30 dyons of 9 kinds (since the rank is 
eight) [17]. 

Returning to the G{2) group, let us introduce 
the numbers vi 2 characterizing the holonomy 
(108): 

^Si = Ui = -M3 - \/3m8, (112) 

= V2 = 3/^3 + VSfls, 

^So = (2 - 21/1 - z/2). 

We imply that So, 1,2 > 0; if ^3,8 are such that 
this condition is violated, one can always perform 
a time-dependent gauge transformation which 
shifts )U3,8 to the 'fundamental' domain where the 
right hand sides of the above equations are non- 
negative. 

We now build the semiclassical vacuum as an 
ensemble of 3 kinds of dyons (and anti-dyons) sat- 
isfying the neutrality condition Ki = 2K2 = 2Kq. 
This is done in the same way and in the same ap- 
proximation as for the SU{N) gauge theory, see 
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or vi = V2 
1 



Ms = 7=- 



(114) 



sections 8-11. The main conclusion there is that 
the free energy of the ensemble is proportional to 
the geometrical mean of the individual dyon ac- 
tions. In the G(2) case the action of the dyon 
of the 1st kind needs to be squared as it enters 
twice more often than the other two, to satisfy the 
neutrality condition. It results in the nonpertur- 
bative free energy induced by dyons (c/. Eq. (69)) 

F = {vlv^(2 ~ 2ui - V2)Y ■ (113) 

The minimum of this expression in 1/1^2 or, equiv- 
alently, in fi^ g, is achieved when the actions for 
all three kinds of dyons are equal, Sq = Si = S2, 
2vi- V2 = \, or 
1 

This is the value of the holonomy that is dynam- 
ically preferred by the ensemble of dyons. 

The final step is to substitute Eq. (114) into 
the Polyakov line (108), and we find 

7r 47r Stt 

TrL = l + 2cos — h2cos( — )-|-2cos( — ) = (!) 
6 6 6 

We conclude that the confining (zero) value of 
the trace of the Polyakov line follows from the 
minimization of the free energy induced by dyons, 
although there are no a priori symmetry reasons 
why this holonomy is privileged. 

17.2. Deconfinement phase transition in 

G(2) 

As temperature increases, the perturbative free 
energy, also a function of the holonomy, starts to 
play an increasingly important role. The pertur- 
bative potential for any gauge group is given by 
the equation 
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where A„ are the 'charged' gluon masses in the 
background of a constant A/^. More precisely, A„ 
are the eigenvalues of in the adjoint represen- 
tation. In the G(2) case there are 12 nonzero 
eigenvalues: 



A,: 



27rT 
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Ms 



± 



V3' 

- [jlj, ± /i8%/3 



±i (M3± 



Ms 
V3 



which should be plugged into Eq. (115). 

The sum of the perturbative (115) and nonper- 
turbative (113) parts of the free energy as func- 
tion of the holonomy /ia^g is plotted in Fig. 16. 
The contour plots are very similar to those of 
Fig. 15 depicting the free energy at different tem- 
peratures for the SU{i) gauge theory. The 12 
symmetrical sectors in Fig. 16 correspond to 12 
different domains for i^i, V2. They can be related 
to each other by time-dependent gauge transfor- 
mations. Going from one sector to another cor- 
responds to changing the roots used to define the 
dyon solutions; physically all of them are equiva- 
lent. 







Figure 16. Contour plots for the potential en- 
ergy as function of the holonomy ^l^ and /ig in the 
(/i3,/j,8) plane at zero {upper left), critical {upper 
right) and high {lower left) temperatures for the 
G(2) gauge group. A zoom at critical tempera- 
ture {lower right) clearly shows the coexistence of 
two minima, typical for the 1st order phase tran- 
sition. 



At temperatures below critical the minimum is 
at a point in the (/ia, /ig) plane that has no obvi- 
ous privilege from the point of view of symmetry. 
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Nevertheless, it is precisely the point where the 
Polyakov loop vanishes, see the previous subsec- 
tion. 

At certain temperature which can be estimated 
through the A scale of the theory a second mini- 
mum starts to develop. At a critical temperature 
it becomes exactly as deep as the first minimum. 
It is the 1st order transition point. Above the 
phase transition the system stays in the second 
minimum which also has no obvious symmetry 
but it now leads to Tr L ^ 0. Eventually at very 
high temperatures the system will move to the 
trivial holonomy at the origin of the contour plots 
in Fig. 16. This picture is in accordance with 
what is observed in recent lattice studies [41]- [44]. 

Qualitatively, exactly the same story happens 
in the SU{'3) gauge theory, see section 16 and 
Fig. 15, with the only unprinciple difference that 
in the SU{3) case the trajectory of the minimum 
in the holonomy plane is along a definite ray, and 
it appears more symmetric to the eye. 

We conclude that the center Zjv symmetry of 
the SU{N) gauge theory is not of fimdamcntal 
importance in the confinement mechanism, as a 
very similar picture of confinement and deconfine- 
ment emerges in the ccntcrlcss G{2) gauge the- 
ory. Dyons, however, are in a position to explain 
both gauge theories as in both cases they dynam- 
ically drive the system to the confining value of 
the holonomy. 

18. CONCLUSIONS 

In these lectures, we have reviewed some es- 
sential tools used to study Yang Mills theories at 
zero, as well as at nonzero temperatures, most 
of which are related to the underlying topology 
of gauge theories. In particular, we have briefly 
outlined 

• instantons in Quantum Mechanics and 
Yang-Mills (YM) theory 

• the Hamiltonian or Schrodinger formulation 
of the YM theory 

• periodic generalization of instantons to 
nonzero temperatures 



• topological classification of YM classical so- 
lutions 

• Bogomolny-Prasad-Sommerfield 
monopoles, also called dyons in the pure 
YM theory 

• the construction of dyons for arbitrary 
gauge groups 

• instantons with a nontrivial holonomy (the 
KvBLL calorons). 

We have then moved to an uncharted terri- 
tory in an attempt to build a semiclassical pic- 
ture based on the ensemble of dyons, that would 
give an explanation of the puzzling phenomena 
associated with the confinement of quarks. In 
particular, we have discussed 

• integration measure over the collective co- 
ordinates of many dyons 

• reduction of the dyons' partition function 
to that of an exactly solvable 3d quantum 
field theory 

• why the 'confining' holonomy is dynami- 
cally preferred by the ensemble of dyons 

• the calculation of the 'electric' string ten- 
sions for any nonzero iV-ality representa- 
tion from the correlation function of the 
Polyakov lines 

• the 'magnetic' string tension (equal to the 
'electric' one despite the lack of apparent 
Euclidean symmetry in the 3d formalism 
used) from the area law for large Wilson 
loops 

• the cancelation of perturbative gluons in 
the Stefan-Boltzmann law in the confining 
phase 

• an estimate of the critical deconfincment 
temperature in units of the fundamental 
string tension 

• confinement and deconfinement in the cen- 
terless G(2) gauge theory. 
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There are still many points in the presented 
dyon mechanism of confinement, that need better 
understanding if not justification. Therefore at 
the moment it is rather a nonperturbative model 
of strong interactions, than a rigorous theory. 
Nevertheless, the overall picture with its good ac- 
cordance with phenomenological and lattice re- 
sults, even at a quantitative level, looks rather 
satisfactory and hopefully is a step in the right 
direction. 
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